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AN ANALYSIS OF FRACTURE TRACE PATTERNS 
IN AREAS OF FLAT- LYING SEDIMENTARY ROCKS 
FOR THE DETECTION OF BURIED GEOLOGIC STRUCTURE 

by 

Melvin H. Podwysocki 
ABSTRACT 

Two study areas in a cratonic platform underlain by flat-lying sedimentary rocks 
were analyzed to determine if a quantitative relationship exists between fracture 
trace patterns and their frequency distributions and subsurface structural closures 
which might contain petroleum. Fracture trace lengths and frequency (number 
of fracture traces per unit area) were analyzed by trend surface analysis and 
length frequency distributions also were compared to a standard Gaussian dis- 
tribution. Composite rose diagrams of fracture traces were analyzed using a 
multivariate analysis method which grouped or ’’clustered” the rose diagrams 
and their respective areas on the basis of the behavior of the rays of the rose 
diagram. 

Analysts indicates that the lengths of fracture traces are log-normally distributed 
according to the mapping technique used in this paper. Deviations from log- 
normality may be associated with both reef (passive) structures whose "closure” 
is caused by differential compaction of sediments over the reefs and with base- 
ment uplift (active) anticlinal-structures. The primary control of fracture trace 
frequency and log-mean lengths is associated with variations in surficial lithology. 
This variation may be extracted using trend surfaces and the residuals may be 
analyzed. Fracture trace frequency appeared higher on the flanks of active 
structures and lower around passive reef structures. Fracture trace log-mean 
lengths were shorter over several types of structures, perhaps due to 
increased fracturing and subsequent erosion. 

Analysis of rose diagrams using a multivariate technique indicated lithology as 
the primary control for the lower grouping levels. Groupings at higher levels 
indicated that areas overlying active structures may be isolated from their neigh- 
bors by this technique while passive structures showed no differences which 
could be isolated. 
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AN ANALYSIS OF FRACTURE TRACE PATTERNS IN 
AREAS OF FLAT- LYING SEDIMENTARY ROCKS FOR 
THE DETECTION OF BURIED GEOLOGIC STRUCTURE 


INTRODUCTION 

Although linear features on the earth's surface had long been mapped solely on 
topographic and geologic criteria (Hobbs, 1911; Brock, 1957), more of these 
subtle features became apparent as aerial photographic coverage became 
available (Rich, 1928). Since then, airphoto linears have been applied to a 
wide range of topics such as groundwater studies (Lattman and Parizek, 1964; 
Siddiqui and Parizek, 1971), mineralization (Keim, 1962; Kutina, 1969) and 
engineering studies (Parizek and Voight, 1970; Parizek, 1971; Alpay, 1973; 
Benedict and Thompson, 1973). Linears observable on various scales of 
aerial photographs and topographic maps have been utilized extensively in 
regional tectonic studies (Plafker, 1964; Gol'braikh et al. , 1968a; Gold et al. , 
1974). Although several investigators claim that analysis of airphoto linears 
will allow exploration for geologic structures which may bear petroleum 
(Permyakov, 1949, 1954; Blanchet, 1957; Mollard, 1957), few exploration 
techniques have been divulged due to their proprietary nature. This paper will 
discuss some parameters which can be extracted from an airphoto linear study 
for the purposes of exploration for several types of oil and gas traps. 

NOMENCLATURE 


The terms "airphoto linears" or "linears" were used above in order to circum- 
vent the variety of names and non- systematic nomenclature for these topographic 
and photographic expressions. Barton (1933) used the term "topographic lines" 
and Gross (1951) used "topographic linears." Although their maps showed they did 
limit the size of the observed features, no comment was made concerning the 
distribution of their individual lengths. Only recently has attention been paid 
to the scale of observations and size of the features (Nemec, 1970; Gold et al. , 
1974) and until the advent of satellite imagery, there was no convenient format 
for direct observations of the large features. 

Blanchet (1957) categorized his observations on linears observed on aerial 
photographs as "micro- and macrofractures, " dividing hie two categories at 
2. 5 miles (4 km). He claimed, but offered no proof, that microfractures (0. 5 - 
2. 5 miles (0. 8 - 4 km)) in length are intrinsic to the sediments themselves 
whereas macrofractures (greater than 2. 5 miles (4 km)) are related to deep 
seated basement features. Because similar orientations prevailed in different 
parts of the world, he claimed that the fractures were related to a worldwide 
tectonic pattern. 
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Mollard (1957) used the term "lineament" to classify aerial photographic linears. 
His classification allowed the use of both continuous and discontinuous features 
ranging from 0. 2 - 5 miles (0. 3 - 8 km) in length. He too considered them 
related to global tectonics. 

Gol'braikh et al. (1968 a,b) use the term "megajoint, " which they adopted 
because of its relationship in hierarchy to other scales of jointing (i. e. micro- 
and macro jointing) and to the analytical techniques which could be applied 
regardless of scale. The term megajoint is based on the scale of maps or 
aerial photographs used and the minimum length (1 cm) which they believe can 
be precisely measured to determine the bearing of a megajoint. Their published 
works indicate a range of 1 to 6 km with a peak around 3 km (Mirkin, 1973, 
pers. comm. ). Unfortunately, this scheme is dependent upon the scale of maps 
or photographs used. Other Russian terms used to describe the same phenomena 
are "lineamental jointing, " "rectilinear elements of topography and stream 
networks" (Gol'braikh et al. , 1968a) and "lineaments" (Shul'ts, 1969). 

Lattman (1958) subdivides airphoto linears into "lineaments" and fracture 
traces," based on their length. He defines fracture traces as naturally 
occurring linear features observed on aerial photographs as alignments of 
stream segments, topographic features and soil and vegetational tonals which 
are expressed continuously for less than one mile in length. He relates them 
either to small faults or zones of joint concentration which are usually vertical 
or nearly vertical in cross-section (Lattman and Matzke, 1961). Excluded from 
the definition are bedding planes, compositional layering, and foliations. 
Lineaments are defined as consisting of the same morphological landscape 
elements as fracture traces, except that they are expressed diseontinuously in 
the landscape and are greater than one mile (1. 6 km) and up to several tens or 
hundreds of miles (km) in length. They may consist of zones of increased 
fracture trace concentrations, transgressing structural, temporal and physio- 
graphic provinces and because of their great lengths, they are thought to be 
recurrent effects associated with basement faults or zones of tectonic adjust- 
ment between major crustal blocks (Wise, 1968; Gold et al. , 1973, 1974). A 
plot of aerial photographic linears combining both of Lattman' s categories 
indicates a bimodal distribution, with a minimum occurring at about the one 
mile length (Lattman, 1969, pers. comm.). Mirkin (1973, pers. comm.) 
indicates a similar bimodal distribution with his break occurring at the 3-4 km 
interval. 

The present study will use the terminology of Lattman (1958) and will examine 
whether his definitions agree with observations made during this study. 
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MEASUREMENT PARAMETERS 


Griffiths (196 7) characterizes the measurable properties of an object by the 
following mathematical equation; 

P = f (material, size, shape, orientation, packing) (1) 

Size (length) and orientation (bearing) are the most readily measured properties 
of fracture traces. Shape can be variously defined. Griffiths (1967) characterizes 
the shape of quartz grains or pebbles as the ratio of their long, intermediate 
and short axes. In this sense, the ratio of fracture trace width to its length 
might be a measurable parameter. However, measurement of fracture trace 
widths is a highly subjective study because of possible erosional and seasonal 
vegetal enhancement, and until more is known of their character with depth, 
no consistent classification can be attempted. In addition, since fracture 
traces are defined as lines, their width can be defined as infinitely small and 
unmeasurable. A radius of curvature can also be defined as a shape parameter, 
however, the scarcity of these features would preclude their use as a commonly 
measured and quantified parameter (Gol'braikh et al. , 1968a). The possible 
significance of these curved and arcuate features has often been overlooked 
(Podwysocki and Gold, 1974); they may represent the surface expression of 
periclinal structures, listric faults and intrusive bodies. 

The two remaining factors which can be studied are materials and packing. In 
this study, materials will refer to surficial geologic materials (formations) 
present in the mapping area. Packing (density or number of fracture traces per 
unit area) will be one of the parameters calculated as a result of this investigation. 

STUDY AREAS 


Two study areas were chosen representing different types of "structural traps" 
for the accumulation of petroleum. Both are located in the relatively stable 
cratonic platform areas of the central USA. A study area in south-central 
Kansas was chosen because it was regarded as typical of vertical uplift 
controlled by basement faulting. The other area was located in west Texas and 
is underlain by a series of reef structures with overlying sediments draping 
over them (differential compaction). No basement tectonic control is evident 
in the latter area. 

The Kansas study area, covering approximately 150 square miles (270 sq. km), 
occupies the southern portion of Pratt and the northern part of Barber Counties 
(Figure 1). It overlies a portion of the southward plunging nose of the Pratt 
Anticline, a southerly extension of the Central Kansas Uplift (Merriam, 1963). 
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KANSAS STUDY AREA 
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Figure 1. Schematic Geologic Map Showing Anticlinal Structures and If 
Productive, the Name of the Associated Oil Fields. The 
Outlines of the Fields Are Based on the Lowest Structure 
Contour Which Indicates Closure, and Structural Closure is 
the Difference Between the Structural Crest and the Lowest 
Contour of Structural Closure. 
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Although deformation occurred as early as Cambrian time (Williams, 1968), 
the major pulse is Mid- Pennsylvanian (Merriam, 1963), and produced an 
unconformity between pre-Mississippian and late Pennsylvanian rocks. 
Structural and structural-stratigraphic traps suitable for the entrapment of 
petroleum were created by "crenulations" of 2-3 km diameter on the Pratt 
Anticline. A northeast trending fault underlines the Coats Oil Field, cutting 
the Precambrian basement (Cole, 1962). No documentation exists for this 
fault in higher stratigraphic horizons (Williams, 1968). Figure 1 contains 
a schematic representation of the oil fields in the area and the amount of 
structural closure as determined by structure contours on top of the Late 
Pennsylvanian Lansing Group (Williams, 1968). According to cross-sections 
by Curtis (1956), minor reactivation of some of the structures may have 
occurred as late as Permian time. Average depth to the top of the producing 
horizons is approximately 3500 feet (1065 m) below the surface and depth to 
basement averages about 5200 feet (1585 m). 

Figure 1 also contains a surface geologic map. Glacial outwash gravels, sands, 
silts and some clays of the Pleistocene Kansan and Illinoian Stages predominate. 
The Illinoian materials are found on the upland surfaces in the northern portions 
of the study areas whereas Kansan materials are usually found in the southern 
part and the major stream valleys of the central portion (Layton and Berry, 
1973). Thickness of these deposits reaches a maximum of 200 feet (61 m) in 
the northern part of the study area and gradually tapers to a zero edge where 
the Permian rocks of the Whitehorse Group crop out in the southern extremity 
of the study area* The latter consist of reddish-brown siltstones, shales and 
sandstones with lesser amounts of gypsum, salt, anhydrite and limestone 
(Layton and Berry, 1973). 

The Texas study area, covering approximately 180 square miles (324 sq. km), 
is situated in the northwestern portion of Nolan and southwestern part of Fisher 
Counties (Figure 2). It lies on the eastern shelf of the Midland Basin, the site 
of the Pennsylvanian and pre- Pennsylvanian Concho Arch and Platform (Hope, 
1956). Two major unconformities exist with a hiatus from Late Ordovician 
through the Mississippian and another from Triassic through the Cretaceous 
ages (Hope, 1956; Shamburger, 1967). During Pennsylvanian time the area 
was the site of extensive reef-building, caused by repetitive advances and 
retreats of the seas across this shallow platform area (Van Siclen, 1958), 
Subsequent deposition commonly covered the reefs with fine-grained clastic 
sediments, eventually draping over them, due to differential compaction, to 
create ''structural highs" (Conselman, 1959), In addition, stratigraphic traps 
associated with the updip pinchout of fore-reef detritus are common. No 
documentation exists for faulting in the study area (Hope, 1956), Depth to 
"Canyon Reef" production horizons averages about 6000 feet (1830 m). 
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TEXAS STUDY AREA 
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Figure 2. Schematic Geologic Map Showing Known Reef "Structures" and 

the Name of the Associated Productive Oil Fields. The Outlines 
of the Fields Are Based on the Lowest Structural Contour Which 
Indicates Closure, and Structural Closure is the Difference 
Between the Structural Crest and the Lowest Contour of 
Structural Closure. 
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Figure 2 contains a schematic representation of reef production with a minimum 
"structural closure" indicated for each reef and excludes stratigraphic traps 
such as the "Canyon Sands, " 

The surface geology is also portrayed in Figure 2. The Permian Whitehorse 
Group consists of sandstone, siltstone and shale redbeds with some inter- 
spersed gypsum beds. It crops out in the northern and extreme eastern 
portion of the study area. The Triassic Dockum Group crops out sporadically 
in the northern and eastern parts of the study area because of its cover by the 
Cretaceous and Tertiary units and its erosion during a later hiatus (Conselman, 
1959; Shamburger, 196 7). This unit consists mainly of red and tan conglomerates, 
sandstones and shales. Because of its small extent and similarity in lithology 
to the Permian, the two units have been grouped together. Dips on the Permo- 
Triassic and older subsurface units is about 0. 5 - 1 degree to the west. The 
Cretaceous Trinity Group occupies the extreme southern part of the study area 
and consists of medium to coarse-grained quartz sands up to 80 feet thick which 
vary in color (Shamburger, 196 7). Directly above is the Fredricksburg Group, 
consisting of thin to thick bedded arenaceous and fossiliferous limestones. 
Maximum thickness approaches 200 feet (61 m) in the Edwards Plateau directly 
to the south of the area, but it is considerably thinner locally. Karst, features 
such as broad, shallow, poorly defined sinkholes are also found. Although 
these limestones underlie most of the central and western portion of the area, 
they are masked by a relatively thin cover of Tertiary Ogallala deposits. The 
Cretaceous units have been consolidated into one map unit because of the small 
lateral extent of the Trinity Group. Regional dips on these units usually do not 
exceed 1 degree to the southeast. The Pliocene Ogallala Formation consists 
of caliche, sands, gravels and some light colored clays; it forms a thin mantle 
over the central and western parts of the field area and, where exposed in 
cross-section in limestone quarries, it does not exceed 8 feet in thickness. 

MAPPING METHOD 

U. S. Department of Agriculture aerial photographs at a scale of 1:20, 000 taken 
in the early 1960’s were used as a basis for mapping the fracture traces. A 
pocket stereoscope was used in areas of moderate relief (up to 150 feet (46 m)), 
and for low relief areas (5-20 feet (1. 5 - 6 m)), individual photographs were 
viewed at low oblique angles while the photos were rotated to view all possible 
"look directions. " Mapping was done in flightlines, spending about 1/2 hour 
per stereo pair. Trainer (196 7) showed that 84-89% of the fracture traces 
could be found in the first 20 minutes of observation. 

As a check to determine if this operator was consistent in the selection of 
fracture traces, parts of the sidelap between adjacent flightlines were mapped 
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and compared. A minimum of 83% of fracture traces were mapped consistently 
between several pairs of flightlines. 

As a test of variations in the recognition of fracture traces, several experienced 
operators were compared to determine if the same general trends were mapped 
amongst the operators. Four sets of airphoto stereo pairs representing 
different types of topography in the study were mapped by two additional 
operators. Freidman Two-Way Analysis of Variance (Siegal, 1956) indicated 
that each operator mapped a different number of fracture traces on the four 
examples, based on a 0. 05 level of rejection. However, the relative ranking 
by the operators of fracture trace direction indicated that in three out of four 
cases, there was no reason to reject the hypothesis that the operators were 
choosing the same directions. Thus, even though absolute numbers of fracture 
traces varied between operators, the same patterns of orientation and the 
same relative magnitude remained when the data were plotted in rose diagram 
plots. GoFbraikh et al. (1968a) achieves the same end by converting the 
absolute number or length of megajoints to percent rose diagrams in order to 
eliminate variation due to different operators and to more clearly discern the 
signal pattern. 

Fracture traces were mapped directly onto aerial photographs by marking 
their endpoints with a soft colored pencil. In order to minimize planimetric 
errors, the fractures were mapped only within a three inch radius of the 
photograph centerpoint. These data were then transferred using a Saltzman 
projector to standard U. S. Geological Survey 1:24, 000 scale topographic maps 
which were used as a base map. Figures 3 and 4 represent the fracture maps 
of the two areas. The grid on the left and top margins will be discussed later. 

Cultural features such as pipelines and fencerows were usually readily distinguish- 
able on aerial photographs. Subsequent field examinations verified and eliminated 
these features. Difficulty was encountered in differentiating some cultural 
features from fracture traces, notably relict plow patterns. This manifested 
itself in two fashions: 1) Plow patterns which paralleled some fracture traces 
would most likely cause the operator to overlook these fractures. This would 
eliminate north-south and east-west oriented fractures in the Kansas area. In 
west Texas, due to the orientation of the cultural pattern, those fractures 
oriented within several degrees of N12W and N78E could be easily overlooked. 
Conversely, old plowing practices did not heed the ’’lay of the land, M and 
plowing was done normal to local slope. Those plow furrows normal to the 
slope would enhance and concentrate runoff in this direction, creating a series 
of parallel first and second order stream channels. Contour plowing practices 
alleviated this problem, however, they may have additionally obscured some of 
the original fracture pattern. Figures 5 and 6 are obvious examples of some 
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FRACTURE TRACE MAP 
S. PRATT AND N BARBER COS. , KS. 



Figure 3. Fracture Trace Map of Kansas Study Area. The Upper and Left 
Margins. Correspond to the Same Margins of Figure 1 and All 
Later Maps Using the Same Base. 
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FRACTURE TRACE MAP 
NW NOLAN & SW FISHER COUNTIES, TEXAS 



Figure 4, Fracture Trace Map of Texas Study Area. The Upper and Left 
Margins Correspond to the Same Margins of Figure 2 and All 
Later Maps Using the Same Base. 
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REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 



Figure 5. Vertical Aerial Photograph of a Portion of Pratt County, Kansas, 
Taken in 1950. Note the Finishing Passes in the Plow Pattern (A), 
the Fracture Trace (B) and Areal Extent of Exposed Carbonate-Rich 
”B" Soil Horizon at C. Compare with Figure 6. 



Figure 6. Vertical Aerial Photograph of a Portion of Pratt County, Kansas, 
Taken in 1963. Finishing Plow Patterns (A) Might be Mistaken for 
Fracture Traces. Fracture Trace (B in Figure 5) Has Been 
Obliterated by Land Contouring. Poor Agricultural Practices Have 
Caused Erosion and Exposed More Carbonate- Rich "B" Horizon (C). 


11 


of these phenomena. Many others exist where a decision concerning their 
origin is more difficult. 2) Plow patterns with their characteristic finishing 
passes through the field diagonals create linear patterns which later show 
through as relict patterns through a newer plowing pattern. Because most of 
these lines pass through field or section corners, they were regarded with 
suspicion and their significance downgraded. Gol'braikh et al. (1968a) noted 
similar problems in the USSR. 

For the purpose of this study and to eliminate some of the subjectivity of the 
mapping, only continuous features were mapped as individual fracture traces. 
Thus, if a linear feature of 4 cm length on an aerial photograph appeared to 
have a break in its length, dividing it into two individual fractures, it would 
be mapped as such. 

DATA HANDLING 

Due to the large amount of information obtained, a computer-based data 
handling system was devised. A cartesian coordinate system was established 
with its origin in the upper left corner of each of the map areas. The X axis 
was chosen as latitudinal and positive to the right and the Y axis meridional 
and positive downward. The beginning and end points of each fracture trace 
could now be referenced with respect to this system which is illustrated in 
Figures 3 and 4. The map data were digitized onto standard 80 column 
Hollerith computer cards and preliminary treatment performed by a FORTRAN 
IV program TRANSFORM. Program listings and additional detail are described 
in Podwysocki (1974). The punched card output of this program contained the 
beginning, end and midpoints of each fracture trace as well as its length in 
millimeters on the map, azimuth and several other parameters which were 
then used in additional computer programs. Subsequent programs utilized the 
established cartesian base, dividing the map area into various grid cell sizes, 
and summarized the data in several fashions. 

These programs were designed so that not only could data be summarized 
within a grid cell specified by the user, but the increment by which this grid 
cell was moved across the map could be specified. Thus, 1) the whole map 
could be treated as a single grid cell and all information would be summarized 
wdthin that one cell, 2) the map could be subdivided into a series of smaller 
cells with the summaries taking place in those individual cells or 3) the map 
could be subdivided as in 2 above and the summary cell size could be incremented 
at a value less than the grid cell size, creating a "running average" or smoothing 
effect (see Podwysocki, 1974). Gol’braikh et al. (1968a) used Hie latter tech- 
nique to look for changes in the number of megajoints and their orientation which 
might be associated with the presence of structural complications (i. e. 
structural closures, faults). 
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ANALYSIS OF FRACTURE TRACE LENGTHS 


Treating the whole map of each area as a single grid cell, and classifying the 
fracture traces into 0. 05 mile (0. 08 km) class intervals, produced the results 
shown in Figures 7 and 8. VECLEN, the computer program for this classifi- 
cation, which is described in Appendix A, summarizes a fracture trace by its 
length if its midpoint falls within a grid cell. In both study areas the distri- 
butions of fracture trace lengths are highly skewed towards the shorter lengths. 
Gol'braikh and Mirkin (1973, pers. comm. ) showed similar results for their 
studies of the Vilyuisk Syneclise and the Preverkhoyansk Downwarp. Although 
no conscientious effort was made by the operator to discriminate against 
linear features greater than one mile (1. 6 km) in length, it should be noted 
that all but a few fracture traces mapped were less than the maximum defined 
length of one mile as defined by Lattman (1958). 

Because of the marked similarity between the observed distribution of fracture 
trace lengths and plots of sediment grain size distribution from sieve analysis, 
a variation of Krumbein's Phi scale transformation (1938) was applied to the 
data as follows : 


z = log 2 x + 6 (2) 

where x is the original length of the fracture trace in miles, z is the transformed 
value of the fracture trace length and 6 is a constant added to each value so that 
all resultant values in this work would be positive. Repeated analysis using the 
same techniques listed above produced the results illustrated in Figures 9 and 
10. The histograms look like Gaussian distributions, however, the summary 
statistics in the figures do not bear this out. The following discussion of 
fracture trace lengths will utilize the transformed data. 

It was thought that mixing of geologically different populations might cause the 
deviations from log-normality in the transformed data. The study areas were 
divided into quarters and each analyzed independently. Results indicated that 
only some areas showed normal distributions. It was noted that the log-mean 
fracture length was different for each of the 4 quarters of each of the two study 
areas. 

To isolate those areas which were anomalous, the study areas were a gain 
quartered, producing a i/16th unit of the total map area and the analysis per- 
formed on each unit. In addition, the summary unit cell was incremented by 
1/2 cell intervals in both the X and Y directions, creating a running average 
as described earlier. The summaries produced cells which were approximately 
3. 5 by 3. 9 miles (5. 6 by 6. 2 km) in the west Texas area and 3. 6 by 2. 9 miles 
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Figure 7. Plot of the Distribution of Fracture Trace Frequency versus Linear Length for the 
Kansas Study Area and the Test for Normality of the Distribution 
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Figure 8. Plot of the Distribution of Fracture Trace Frequency versus Linear Length for the 

Texas Study Area and the Test for Normality of the Distribution 
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Figure 9. Plot of Fracture Trace Frequency versus Log- Length for the Kansas Study Area 
and the Test for Log-Normality of the Distribution. Numbers Within Parentheses 
Represent Values When Distribution Tails Were Folded So That Expected 
Frequency > 0. 95. 
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Figure 10. Plot of Fracture Trace Frequency versus Log- Length for the Texas Study Area 

and the Test for Log-Normality of the Distribution. Numbers Within Parentheses 
Represent Values When Distribution Tails Were Folded So That Expected 
Frequency > 0. 95. 


Reproducibility of this 

ORIGINAL PAGE IS POOR 



(5. 8 by 4. 6 km) in the Kansas area for a total of 64 cells (8 by 8 in each area). 
Summary statistics such as the log-mean fracture trace length, standard 
deviation, skewness (\Zj3j ) and kurtosis ( jS 2 ) of each cell’s frequency 
distribution as well as the number of fracture traces for each unit cell 
were produced for use in additional analyses. 

The summary statistics produced in the above mentioned compilations of the 
data were analyzed using linear regression analysis. Due to the paucity of 
fracture traces in the southernmost tier of cells (less than 5 in each) in the 
Texas study area, these cells were eliminated from the analysis. 

The significance test for correlations between the statistical moments for the 
Kansas data (Table 1) indicates a significant correlation for 1) log-mean 
fracture trace length and skewness, 2) number of fracture traces per unit cell 
and standard deviation and 3) skewness versus kurtosis. Figure 11 represents 
the plot of the standard deviation versus number of fracture traces per unit 
cell. The plot indicates low standard deviations associated with cells containing 
few fracture traces (lower left part of diagram). Because the reliability of the 
statistical moments for such small sample sizes is highly questionable, the 
offending samples (all cells containing less than 45 samples), which occurred 
along the eastern and southern margins of the map area, and were due to 
incomplete mapping coverage, were eliminated from consideration in further 
tests. Repeated regression analysis on the data exclusive of the mentioned 
marginal cells indicated no significant correlation between two of the three 
previously determined associations. However, it should be noted that a 
significant correlation did remain between the skewness and kurtosis measures; 
Figure 12, based on the original analysis of 64 Samples, serves to illustrate 
the results. A small group of samples located near the right margin of the 
plot contains kurtosis values which are highly leptokurtic * (8-10). These cells 
contain several very long fractures (greater than the accepted length for a 
fracture trace) that were inadvertently included, and will be discussed later 
in the log-normality analysis. A removal of these four anomalous cells and 
repeated regression analysis indicated no significant correlation between the two 
moments. Removal of these correlations, or attributing them to some sampling 
inconsistencies, indicates the samples are homogeneous, that is, several 
discrete and very distinct populations do not exist in the data. 


* More peaked lhan normal 
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Table 1 


Results of Linear Regression Analysis 
On Log-Mean Fracture Trace Moments 


Kansas Data - 64 Samples 



Log-Mean 

Length 

Standard 

Deviation 

Skewness 

Kurtosis 

Standard 

Deviation 

NS 




Skewness 

s* 

NS 



Kurtosis 

NS 

NS 

S** 


No. of Fracture 
Traces per 
Unit Cell 

NS 

S** 

NS 

NS 



Table 2 




Results of Linear Regression Analysis 
On Log-Mean Fracture Trace Moments 



Texas Data - 56 Samples 




Log-Mean 

Length 

Standard 

Deviation 

Skewness 

Kurtosis 

Standard 

Deviation 

g** 




Skewness 

g** 

NS 



Kurtosis 

NS 

S* 

NS 


No. of Fracture 
Traces per 
Unit Cell 

NS 

NS 

NS 

g** 


NS = non significant 

S* = significant at 0.05 level 

S** - significant at 0.01 level 
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Figure 11. Regression Analysis Plot of Standard Deviation versus the Number of Fracture Traces 
Per Unit Cell for the Kansas Data. Numbers Within the Plot Indicate the Number of 
Data Points Located in that Position. 


REPRODUCIBILITY OP THE 
ORIGINAL PAGE IS POOR 



O Vi 


KANSAS STUDY AREA 


SKEkNESS (ROOT ai I 


K 

C 

a 

i 


( 

f 

2 

I 


- 0*247 


-0.007 


0*332 


0.672 


1*012 


1.3*1 


1*691 


2*031 


2*370 


2*710 


* 


»..**•• *4 ... *6 ..»***«. 


23*399 
22.049 
22*399 
21*899 
21 .399 
20*899 
20*399 

19.899 
19*399 

18.899 
18*399 

17.899 
17.399 
1 0*899 
16*399 
19*699 
15*399 
14.699 
14*399 
13*699 
13*399 
12*699 
12*399 
11*699 
11*399 
10*899 
10*399 

9.699 
9*399 
8*899 
6.359 
7*099 
7*399 
6*699 
6*399 
5*699 
5*399 
4.899 
4*399 
3*699 
3.399 
2*899 
2*399 
1 .899 
1 .399 



-0.007 


0.672 


1*351 


2*031 


2*710 


23.399 

22.899 
22*399 
21*899 

21.399 

20*699 

20.399 
1 9* 899 
1 9* 399 
I 8*899 
l 8. 399 
17*899 

I 7. 399 

16. 899 
16*399 

15.899 

15.399 

14. 899 

14.399 

1 3.899 

13.399 
12*699 
12*399 

I I * 899 
1 1*399 
1 0*899 
10*399 

9.899 
9*399 
8.099 
0.399 
7* B99 

7.399 
6* 899 

6.399 
5. 899 
5. 399 

4.899 
4*399 
3*899 
3*399 
2*899 
2*399 
1*899 
1* 399 


O & 

IS 

2 § 

£3 £4 


w o 

tl 


Figure 12. Regression Analysis Plot of Skewness versus Kurtosis for the Kansas Data. 



The regression analyses on the Texas data are given in Table 2. Significant 
correlations exist between: 

1) log-mean fracture trace length and standard deviation 
(Figure 13), which indicates increasing standard deviation 
with increasing mean fracture trace length; 

2) log-mean fracture trace length and skewness (Figure 14), 
which illustrates an increasing positive skewness (mode 
displaced towards smaller values with respect to the 
mean of the distribution) with increasing fracture trace 
length; 

3) standard deviation and kurtosis, (Figure 15); and 

4) kurtosis and the number of fracture traces per unit cell 
(Figure 16). 

Figures 13-15 can be interpreted together to indicate one of two possible 
causes. If the assumption is made that the samples were taken from a single 
homogeneous population, then the sampling technique indicates a bias. Con- 
versely, the population may not be homogeneous, and the tests may indicate 
the sampling of two or more discrete and distinct populations of fracture 
traces. The second of the two hypotheses will be proven and more clearly 
illustrated by the use of trend surface analysis which will be discussed later. 

The last significant correlation occurs between kurtosis and the number of 
fracture traces per unit cell (Figure 16). These high values are associated 
with large sample populations and are anomalous, perhaps suggesting some 
mixing of several populations of fracture traces. 

Tests were performed using the Chi Square, skewness and kurtosis criteria 
(Griffiths and Ondrick, 1968), comparing the observed against a hypothetical 
Gaussian distribution. Deviations of each of the criteria were ranked, 
assigning values to those populations which significantly differed from 
normality at the 0. 05 and 0. 01 levels. Rankings were assigned as illustrated 
in Table 3. 

If a criterion value was non- significant, it was assigned a zero value. The 
rankings of the three criteria for each cell were then summed to create an 
index value characteristic of the population distribution in each cell. High 
ranking values indicate strong deviations from log-normality as illustrated 
in Figures 17 and 18. 
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Figure 13. Regression Analysis Plot of Log-Mean Fracture Trace Length versus Standard 

Deviation for the Texas Data. 
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Regression Analysis Plot of Standard Deviation versus Kurtosis for the Texas Data, 
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Cell for the Texas Data. 
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Figure 17. Results of Test for the Distribution of Fracture Trace Lengths to 
Log-Normality. Rank-values Are Associated With the Centerpoint 
of Each Grid Cell. Refer to Table 3 and Text for Key to Rank 
Values. 
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Figure 18. Results of Test for Distribution of Fracture Trace Lengths to 
Log-Normality. Rank-values Are Associated with the Center- 
point of Each Grid Cell. Refer to Table 3 and Text for Key to 
Rank Values. 
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Table 3 


Rankings for Deviations from Log-Normality 


Criterion 

Level of 
Significance 

Rank 

Chi Square 

0. 05 

2 


0. 01 

4 

Skewness 

0. 05 

1 


0. 01 

2 

Kurtosis 

0. 05 

1 


0. 01 

2 


In many cases, significant deviations from log-normality occur over known 
structures. Analysis of fracture trace log-lengths in the Kansas study area 
indicates a consistent positive skewing (mode displaced towards smaller values 
with respect to the mean) in cells which rank four or higher; their respective 
kurtosis values are leptokurtic. Several factors may account for these varia- 
tions. First, structural control may exist, possibly causing development of 
shorter fractures over structures due to enhancement of surface factors such 
as erosion along the fractures. Secondly, control may be due to changes in 
lithology. Analysis of fracture trace lengths does indicate a lithologic control 
and will be discussed shortly. Thus, mixing of two surface rock types within 
a grid cell may cause this type of discrepancy. However, it should be pointed 
out that similar mixing also takes place in the two Pleistocene aged formations 
in the eastern part of the study area, and these types of deviations do not exist 
in this area. Very high deviations in the southern portion of the Kansas area 
(ranked 7 and 8) may be due to the proximity to Permian outcrops and/or the 
influence of two fractures greater than one mile (1. 6 km) in length that were 
inadvertently mapped (see Figure 7). Because these features were larger by 
a factor of two over all other fracture traces in the area, they cause highly 
significant deviations from log-normality. These same fractures were re- 
sponsible for the high correlation between the skewness and kurtosis in the 
regression analysis plot (see Figure 12). Thirdly, biases due to operator 
fatigue or cultural land practices may occur. These hopefully were minimized 
with field checking, rest periods during mapping and cross checking with 
photographic coverage of earlier dates to eliminate these possible errors. 

Significant deviations from log-normality (ranked four or higher) also occur 
over some of the known reefs in the west Texas study area. Skewness and 
kurtosis behave similarly to the anomalies in the Kansas area. The same 
three arguments stated in the previous paragraph may be employed. Cultural 
effects have been minimized by reference to earlier photographic coverage. 
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Because lithology shows little control in the northern part of the area where 
cells transgress lithologic boundaries, it is probably not a controlling factor 
in the anomalous eastern portion of the study area. The high value (a rank of 
6 in row 4, column 3) in the central portion of the map area is due to the 
presence of a lineament. The fracture trace length distribution for this area 
is unlike those over the reefs; it is skewed positive and is nearly normal in its 
kurtosis. In some cases, the anomalous ranks do not directly overlie the 
structure, but lie on its flanks. Harris et al. , (1960); Gol'braikh et al. , 
(1968a) and Saunders (1969) indicate that increased fracture density may occur 
along the flanks of a structure, however, no mention has been made of changes 
in fracture length. 

Further reduction of the grid cell size produced many cells with too small a 
population, and thus reliable statistics were not possible. Analysis of fracture 
trace length distributions in individual 10 degree azimuth classes in each grid 
cell also proved fruitless because of the small number of fracture traces in 
each cell. 

ANALYSIS OF FRACTURE TRACE FREQUENCY 


Trainer and Ellison (1967) define frequency as the number of fracture traces, 
irrespective of their length, which fall within a unit area under consideration. 
Trend surface analysis (O’Leary et al. , 1966) was applied to the fracture trace 
frequency values generated by the VECLEN program for the l/16th unit areas 
discussed above. This technique attempts to fit surfaces which represent 
polynomial equations of increasing order to map data. Increasing polynomial 
order represents increasing complexity of the surface, which thus more closely 
approximates the given data. It can be used in some instances to extract 
different components responsible for variations which may be present in the 
data. In most cases, first through sixth order surfaces were fitted to the data. 
Analysis of variance was applied to the output statistics of this technique to 
determine which surfaces were a significant improvement over their lower 
order neighbors (Krumbein and Graybill, 1965); the probability level used was 
based on P = 0. 005. Only selected surfaces which achieved the prescribed 
level of significance and their residual plots will be discussed. Tables 4 and 5 
summarize the data for each study area. 

Figure 19 illustrates the second order surface for the Kansas data and accounts 
for 82% of the variations. It shows that fracture trace frequency is highest in 
the southeastern part of the map area near the Permian outcrops, and decreases 
northward toward the younger Pleistocene deposits and towards the map 
peripheries, where coverage is incomplete or control is lacking. This suggests 
that lithology may be a controlling factor for one of several reasons. 1) The 
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Table 4 


Analysis of Variance of Trend Surfaces Data for 
Fracture Trace Frequency Kansas Study Area 


Surface Order 

Sum of 
Squares 

Degrees of 
Freedom 

Mean 

Square 

F 

Significance 

Cumulative Percent 
Variation Explained 

Percent Variation 
Improvement for 
Each Surface 

1st 

7238.8 

2 

3619.4 

11.2 

.005-001 

39.8 

39.8 

Dev. from 1st 

10964.8 

34 

322.5 





2nd 

7711.8 

3 

2570.6 

24.5 

<.001 

82.1 

42.3 

Dev. from 2nd 

3253.0 

31 

104.9 





3rd 

781.4 

4 

195.4 

2.1 

.10-.25 

86.4 

4.3 

Dev. from 3rd 

2471.6 

27 

91.5 





4th 

1963.1 

5 

392.6 

17.0 

.01 '.025 

97.2 

10.8 

Dev. from 4th 

508.5 

22 

23.1 





5th 

158.2 

6 

26.3 

1.2 

.25-30 

98.1 

0.9 

Dev. from 5 th 

350.3 

16 

21.9 





6th 

238.1 

7 

34.0 

2.72 

o 

o 

99.4 

1.3 

Dev. from 6th 

112.2 

9 

12.5 






Table 5 

Analysis of Variance of Trend Surface Data for 
Fracture Trace Frequency Texas Study Area 

Surface Order 

Sum of 
Squares 

Degrees of 
Freedom 

Mean 

Square 

F 

Significance 

Cumulative Percent 
Variation Explained 

Percent Variation 
Improvement for 
Each Surface 

1st 

8142.2 

2 

4071.1 

13.0 

<.001 

37.7 

37,7 

Dev. from 1st 

13446.7 

43 

312.7 





2nd 

6432.7 

3 

2144.2 

12.2 

<001 

67,5 

29.8 

Dev. from 2nd 

7014.0 

40 

175,4 





3rd 

1041.4 

4 

260.4 

1.6 

.10-25 

72.3 

4.8 

Dev. from 3rd 

5972.6 

36 

165.9 





4th 

3443.6 

5 

688.7 

8.4 

<.001 

88.3 

16.0 

Dev. from 4 th 

2529.0 

31 

8 L58 





5th 

896.4 

6 

149.4 

2.3 

.05-. 10 

92.4 

4.1 

Dev. from 5th 

1632.6 

25 

65.3 





6th 

908.8 

7 

129.8 

. 3.2 

.01.025 

96,6 

4.2 

Dev, from 6th 

723.8 

18 

40.2 
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Figure 19„ Second Order Trend Surface for Fracture Trace Frequency 
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unconsolidated Pleistocene sediments may have a masking effect, subduing the 
number of fractures propogated to the surface; 2) the younger sediments may 
have been subjected to lower stress levels, fewer periods of deformation and 
a shorter time for the propogation of the fractures; or 3) different rock types 
may have different mechanical properties. Because Pleistocene unconsolidated 
deposits do thicken northward, the first two factors are probably the most 
significant. 

Analysis of the residuals* map (Figure 20) indicates a large positive residual 
(greater than the calculated model) in the southeast, part of the map, underlain 
by outcropping Permian rocks, and may be explained by several factors. 

Harris et al. (1960) noted changes in jointing frequency due to contrasting 
lithologies over the Goose Egg Dome in Wyoming. Not only did they find a 
progressive decrease in frequency from siliceous limestone, calcareous 
quartz sandstone, soft sandstones to ductile shales, but also that fracture 
frequency was inversely proportional to strata thickness. In his study on joints 
in the Great Scar limestones in England, Doughty (1968) recorded changes 
between differing limestone types of similar age. Huntington (1969) found 
changes in fracture trace frequency due to contrasting lithology and suggests 
that observations be confined to like rock types. DeSitter (1964) recognized 
lithology and strata thickness, amongst others, as controls of rock fracturing 
intensity. Another factor which should be considered is the possible masking 
of the fractures due to the strong contrast in mechanical properties of the 
consolidated Permian deposits as opposed to the unconsolidated Pleistocene 
materials, which could act as a filter, either totally obliterating or subduing 
some fracture traces,. 

Another positive residual is associated with a series of structural closures in 
the vicinity of the town of Coats (Figure 20). Although the anomaly overlies 
two different map units, the mechanical contrast between these two unconsolidated 
Pleistocene deposits should be minimal. Excluding possible operator bias, 
the residual might reflect the subsurface Pennsylvanian structures. Residuals 
along the map peripheries are discounted due to lack of control. Gol'braikh 
et al. (1968a), Saunders (1969) and Dranovskii (1970) have suggested that the 
number of airphoto linears per unit area (frequency) is an indicator of structural 
culmination. Moreover, Dranovskii (1970) further states that in box- like 
uplifts, maximum fracturing occurs on the fold limbs, while in ridge- like 
uplifts it develops on the crest of the structure. 

The second order surface for the Texas data accounts for 67% of the variation 
and is illustrated in Figure 21. It shows fewer fracture traces over the 


* For any given observed data point on the map: residual “ observed - expected value calculated for the 
coordinates of tile observed data points. 
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Figure 20. Map of Second Order Trend Surface Residuals for Fracture 

Trace Frequency 
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Figure 21. Second Order Trend Surface for Fracture Trace Frequency 
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Tertiary deposits and the immediately underlying Cretaceous limestones, 
whereas more fracture traces occur over the Permo-Triassic rocks. In 
conformity with the previously stated conclusions of Harris et al. (1960) and 
other workers, the same reasons may explain the lower frequency over the 
Cretaceous -Tertiary rocks. The Ogallala Formation forms a thin blanket, 
not exceeding 6-8 feet in thickness over the study area. In addition, inspec- 
tion of several quarries in the Fredricksburg Group limestones revealed a 
large population of curved joint surfaces which usually terminated at bedding 
planes. These are non- systematic joints that are not associated with quarrying 
operations. The paucity of vertical systematic joints suggests that most 
stresses may have been taken up and diffused in the non- systematic joints, 
thereby precluding the formation of wide zones of weakness suitable for the 
development of fracture traces. Because the fourth order residual map more 
clearly illustrates the results, a discussion of the second order map residual 
is unnecessary. 

Figure 22 shows the results of the fourth order fit and answers 89% of the 
variation. The model contours tend to parallel the north-south flightlines, 
which suggests an operator bias due to changes in accuity during mapping, 
however, higher frequencies again occur over the Permo-Triassic rocks. 

This inter-flightline variation was the predominating signal in the residual 
plot of the second order surface. It is therefore suggested that mapping of 
fracture traces either be done on a suitable scaled mosaic or that individual 
photographs or pairs should be picked randomly from the total available set 
so that this type of variation might be distributed more evenly. 

Figure 23 contains the residuals map based on the fourth order surface. 

Although some alignment parallel to the north-south flightlines does occur, 
most has been removed by this surface. The large positive anomaly in the 
northern portion of the area is associated with the "saddle" in the trend surface 
(Figure 16) and is anomalous. The strong negative anomalies in the eastern 
portion of the map area appear to be associated with the flanks of three of the 
reef structures. This observation is further enhanced by the fact that they 
occur along several flightlines, thereby indicating a consistency between flight- 
lines after removal of the inter-flightline variation. 

In summary, the predominating portion of the variation in frequency of fracture 
traces is associated with differences in lithology. Lesser amounts of the 
variation are due to operator variability due to changes in perceptibility during 
fracture trace mapping. Another variation which may occur is associated with 
"structural closures. " Basement uplift structures are accentuated by positive 
(high) fracture trace frequencies along their flanks, whereas "passive" structures 
such as reefs, may be associated with negative (low) fracture trace frequencies 
in these study areas. 
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Figure 22. Fourth Order Trend Surface for Fracture Trace Frequency 
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Figure 23. Map of Fourth Order Trend Surface Residuals for Fracture 

Trace Frequency 
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ANALYSIS OF LOG-MEAN FRACTURE TRACE LENGTHS 

Log-mean fracture trace lengths generated for the l/16th unit areas by the 
computer program VECLEN were also analyzed using trend surface analysis. 
Significance of improvement in the information level of each surface was 
tested as described earlier. Tables 6 and 7 summarize the results. Only the 
highest order surface showing the prescribed level of significance (less than 
0. 005) will be discussed. 

Figure 24 illustrates the fifth order surface for the Kansas study area and 
accounts for 80% of the variation. The model shows longer fracture traces in 
the northern part of the study area, becoming progressively shorter towards the 
Permian outcrop area. This may be interpreted as a masking effect of the 
glacial overburden, causing the operator to overlook shorter fracture traces 
due to their less pronounced nature or their complete obliteration by the 
overburden. The model also shows a parallelism between some of the contours 
and geologic formation boundaries, as exemplified by the 4. 1 contour, which 
further reinforces the lithologic control hypothesis. The parallelism of contours 
and their steep gradient in the western part of the area is due to lack of control 
in this area. 

The corresponding residuals map (Figure 25) indicates a broad positive residual 
trending northwest in the central part of the area, and parallels the boundaries 
between the two mapped Pleistocene units. The positive bands in the northeast 
and southwest sectors also may be associated with the formational boundaries. 

A negative residual is present in the west-central portion of the map and 
coincides with the increased fracture trace frequency derived from the second 
order residual (Figure 20). These two factors may be inter- related; increased 
deformation may cause more intense fracturing (higher frequency) and because 
of surfieial processes, greater erosion generates more linear first and second 
order streams, which manifest themselves as fracture traces. 

The fifth order surface for the west Texas study area (Figure 26) indicates 
longer fracture traces over the Cretaceous-Tertiary deposits, with shorter 
fracture traces occurring in the Permo-Triassic rocks. The observed differ- 
ences may be due to masking effects as discussed earlier for the Kansas area, 
however, Trainer and Ellison (196 7) found that longer fractures traces occurred 
in the limestone units of the Shennandoah Valley. They suggested that this 
might be due to solution and coalescence of joint planes and zones of weakness, 
a process which has operated in this area as evidenced by the development of 
karst features. 

Figure 27 illustrates the residuals map associated with the fifth order surface. 
No consistent pattern is found with respect to the reef structures. The dominant 
features include a negative residual trending northwest in the central part of 
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Table 6 


Analysis of Variance of Trend Surface Data for 
Fracture Trace Log-Mean Length Kansas Study Area 


Surface Order 

Sum of 
Squares 

Degrees of 
Freedom 

Mean 

Squares 

F 

Significance 

Cumulative Percent 
Variation Explained 

Percent Variation 
Explained By 
Each Surface 

1st 

.0514 

2 

.0257 

7.79 

<.001 

21.5 

21.5 

Dev. from 1st 

.1872 

56 

.0033 





2nd 

.0696 

3 

.0232 

10.55 

<.001 

50.6 

29.1 

Dev. from 2nd 

.1176 

53 

.0022 





3rd 

.0167 

4 

.0042 

2.00 

.10-. 25 

57.7 

7.1 

Dev. from 3rd 

.1009 

49 

.0021 





4th 

.0230 

5 

.0046 

3.56 

.025-. 05 

67.4 

9.7 

Dev. from 4th 

.0779 

44 

.0018 





5 th 

.0309 

6 

.0052 

4.33 

.005 -.001 

80.3 

12.9 

Dev. from 5th 

.0470 

38 

.0012 





6th 

.0204 

7 

.0029 

3.22 

.01 -.025 

88.9 

8.6 

Dev. from 6th 

.0266 

31 

.0009 





Table 7 


Analysis of Variance of Trend Surface Data for 
Fracture Trace Log-Mean Length Texas Study Area 


Surface Order 

Sum of 
Squares 

D cgrees of 
Freedom 

Mean 

Squares 

F 

c Cumulative Percent 

Significance Variation Explained 

Percent Variation 
Explained By 
Each Surface 

1st 

1.58 

2 

.790 

29.26 

<001 

48.8 

48.8 

Dev, from 1st 

1.66 

61 

.027 




2nd 

.71 

3 

.037 

14.81 

<.001 

70.7 

22.1 

Dev. from 2nd 

.95 

58 

.016 




3rd 

.28 

4 

.070 

5.83 

<.001 

79.2 

8.3 

Dev. from 3rd 

.67 

54 

.012 




4th 

.30 

5 

.060 

7.50 

<001 

88.6 

9.4 

Dev. from 4th 

.37 

49 

.008 





5th 

.14 

6 

.023 

4.60 

< .001 

92.9 

4.3 

Dev. from 5 th 

.23 

43 

.005 




6th 



NOT ANALYZED 



Dev. from 6th 
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KANSAS STUDY AREA 



Figure 24. Fifth Order Trend Surface for Log-Mean Fracture Trace Length 


41 




KANSAS STUDY AREA 



Figure 25. Map of Fifth Order Trend Surface Residuals for Log-Mean 

Fracture Trace Length 
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TEXAS STUDY AREA 



Figure 26. Fifth Order Trend Surface for Log-Mean Fracture Trace Length 
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TEXAS STUDY AREA 



Figure 27. Map of Fifth Order Trend Surface Residuals for Log-Mean 

Fracture Trace Length 
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the map area. Its extension across several flightlines tends to bear out the 
reality of this feature. Examination of an aerial photographic mosaic reveals 
a lineament passing through the area in this direction which is most likely 
related to this anomaly. The same argument previously discussed concerning 
increased structural deformation, which produces a greater number of shorter 
fractures, may be invoked. Again, some variation is noted between flightlines 
and some of the positive anomalies in the south. 

In summary, log-mean fracture trace lengths are predominantly controlled by 
the type of sedimentary material. Variations between flightlines may occur, 
but their effect is not as pronounced as in fracture trace frequency. Fracture 
traces may be shorter over basement uplift structures. This effect may not 
occur over passive structures, such as those formed by a draping of sediments 
over bedrock highs. 

ANALYSIS OF ROSE DIAGRAMS 


Several formats are available for displaying directional data. Three dimensional 
data, such as attitude of joint planes, can be efficiently portrayed on stereo- 
graphic (Wulff) or equal area (Schmidt) nets. Statistical analysis of these data 
are cumbersome, but has been discussed by several works (Chayes, 1949; 
Fisher, 1953; Pincus, 1953). Two dimensional data, such as the strike of 
fracture traces, may be displayed as histograms or as rose diagrams 
(Podwysocki, 1974). Both these formats can be conveniently tested and may be 
generated by computers and will be used in this paper because of their 
suitability for visual comparison. 

Several methods can be utilized to summarize the data for rose diagrams. 

Trainer and Ellison (196 7) use the terms "frequency 1 ' and "density. " Frequency 
as described earlier, refers to the number of fracture traces, irrespective of 
their length, while density refers to the total length of fracture traces. Each 
of these respective techniques has its disadvantages. Summarization using 
frequency eliminates a bias due to length. Thus, a fracture trace of 0. 25 mile 
(0. 4 km) is given as much weight as one 0. 75 mile (1. 2 km) long. However, 
due to the mapping technique employed in this paper, which breaks up fracture 
traces into components based on their continuous exposure, the shorter fracture 
traces are favored. Thus, density was chosen as the analysis criterion in this 
work in order to minimize this bias. In addition, the need for a standard unit 
to facilitate comparison has been noted elsewhere. Gol'braikh et al. (1968) 
suggest conversion of the units into percent values prior to plotting, so that 
the size of all roses will be standardized. 
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Joints were measured in several bedrock exposures in the west Texas study 
area. Their orientation frequencies were compared to the density of rose 
diagrams of fracture traces measured in grid cells approximately three miles 
square surrounding each of these localities. Data in 10 degree azimuth classes 
were anatyzed using AZMAP and ROSE, computer programs written by 
Podwysocki (1974). In order to compare statistically the two dissimilar units 
of measurement, both sets of data were converted to percentages. A total of 
69 systematic joints were measured in Cretaceous limestones of the 
Fredricksburg Group, exposed in a quarry near the central-western edge of 
the map area. Nearly all systematic joints were vertical, eliminating the need 
to use three-dimensional displays and making the measurements suitable for 
comparison with the fracture trace distribution. As discussed earlier, there 
also were many non- systematic joints. Permian exposures in the extreme 
central-eastern part were measured and consisted of a roadcut in a gypsiferous 
sandstone and a railroad cut in a massive sandstone. A total of 59 joints were 
measured and combined from these two adjacent cuts. All joints in these cuts 
were within 5 degrees of vertical. Figure 28 contains a graphical comparison 
of the two sets of patterns. 

Neither set of rose diagrams show a good visual fit; a Chi Square test comparing 
the fracture trace and joint orientations for each locality indicates that the 
patterns were not similar based on a 0. 01 level of rejection. Neither could it 
be influenced that much by population size, because GoPbraikh et al. (1968a) 
indicated that 40 - 50 joint measurements were required to achieve statistical 
reliability. It should be noted that while there is conformity in direction in 
the Permian rocks there is a consistent angular displacement between the two 
patterns in the Cretaceous rocks. The former set may reflect fracture traces 
that are occupied by zones of joints sub-parallel to the direction of the fracture 
trace (Lattman, 1969, pers. comm. ); the latter may represent a displacement 
of the second order joints from the direction of maximum shear stress. This 
phenomenon has been documented by Renner (1969) and may relate to a 
hierarchichal structural framework as postulated by Moody and Hill (1956) and 
discussed by Nemec (1970) and Gold et al. (1973). Because of the dissimilarity 
in the Cretaceous patterns and partial agreement in the Permian patterns, it 
might also be suggested that either the rocks have behaved differently when 
subjected to the same stresses or that the older units were subjected to an 
additional period of stress not experienced by the younger units. 

These results are partly contrary to those of Lattman and Nickelsen (1958), 
Hough (1959), Boyer and McQueen (1964) and Alpay (1973), who generally 
found good agreement between fracture trace and joint directions in their 
investigations in sedimentary rocks dipping less than 5 degrees. Matzke 
(1961), Lattman and Matzke (1961, 1971) and Trainer and Ellison (1967), 
however, reported that fracture traces and joint directions do not totally 
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FRACTURE TRACES 

JOINTS 


Figure 28. Comparison of Joint and Fracture Trace Hose Diagrams, Texas 

Study Area 
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coincide. Their observations were made in more deformed rocks (i. e. the 
Appalachian fold belt). Lattman and Matzke (1961) suggest that joint patterns 
in relatively stable cratonic areas are paralleled by fracture traces whereas, 
local structure in highly deformed materials impress their own local joint 
sets which may deviate from the regional trends. 

Although orientation directions coincide for the Permian rocks, the length of 
the rays (degree of preferred orientation) is greater for the joints. This is 
probably due to the big difference in the size scale of the areas sampled 
(9 square miles (23 sq. km) versus 2 outcrops 1/2 mile apart). 

ANALYSIS OF FRACTURE TRACE PATTERNS 

Pattern recognition of preferi’ed orientation in fracture analysis tends to be 
more difficult due to the large amount of data and its multivariate nature. 

Several approaches have been used to enhance patterns. Haman (1961, 1964) 
isolated and plotted all macrofractures ’(lineaments) and mesofractures 
(fracture traces) which fell within narrow azimuth ranges and used them in a 
qualitative fashion to discern faulting and to locate changes in regime of 
individual tectonic blocks. Maffi and Marchesini (1964) describe the use of 
optical and computer processing techniques to filter and isolate individual 
trends. Gol'braikh et al. (1968a) also isolated regional structures by plotting 
their megajoint densities for narrow azimuth ranges, and showed the applica- 
bility of Permyakov's (1949) "rule of the parallelogram" to determine regional 
trends by analysis of the rose diagram modes. Little has been published on a 
method for the comparison of several rose diagrams. Chudinskii (Mirkin, 

1973, pers. comm. ) suggests that rose diagrams of small subsets of the total 
area should be compared against the grand rose diagram for the whole territory. 

A variation of the Chi Square criterion could then be used to compare the subset 
against the composite rose diagram. Those which proved to vary significantly 
from the composite diagram were zones of "tectonic complications. " Lattman 
(1969, pers. comm. ) suggested a similar technique, but instead of comparing 
a subset against the composite rose diagram, the subset was compared against 
all of its adjacent neighbors. Significant variations between neighboring diagrams 
would then indicate structural complexities. 

The fracture trace data compiled by the TRANSFORM program was processed 
by AZMAP (Podwysocki, 1974), which classified the fracture traces into 
direction categories within each unit cell (1/16 th of the total study area). As 
described previously in the analysis of fracture trace frequency and lengths, 
a 1/2 cell sliding average increment also was used. An azimuth class interval 
of 10 degrees was utilized during the classification. 
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Several classification techniques could be used in AZMAP. The first, entitled 
"Part" analyzed only that portion of the fracture trace length which lies within 
the cell. "Mid" considered the whole fracture trace within the cell if its mid- 
point fell within the cell. A comparison of the two techniques showed that 
there was no significant difference if the results of the two classification tech- 
niques were compared against each other for each of the 49* grid cells of each 
area, using the Chi Square test and a rejection level of 0. 05. 

Punched card output of the summary length of fracture traces per azimuth class 
per grid cell were processed bj' a computer program ROSE (Podwysocki, 1974), 
which produced i’ose diagrams (see Figures 29 and 30). The punched card 
output from AZMAP also was utilized in a multivariate analysis computer 
program CLUS (Rubin and Friedman, 1967). Each rose diagram consisted of 
18 variables or measurements (the sum total length of fracture traces within 
each of the 10 degree azimuth classes). A total of 49 grid cells (objects) were 
generated by AZMAP for each study area and these were treated as 49 samples. 

Multivariate techniques have been shown by Dahlberg and Griffiths (196 7) to be 
an effective method for determining the relationships between objects with 
interacting properties. The Rubin and Friedman program is appropriate for 
determining the relationships between samples because the procedures allow 
classification on the basis of a number of groups determined by the user. A 
determination of the optimum grouping is made on the basis of several computer 
generated criteria for each classification. 

The inverse of the Wilk T s lambda criterion, log (max | T | / | W | ), is used 
as an informal indicator of the best number of groups (Friedman and Rubin, 
1967), where: 

W is the pooled within-group matrix of the cross products of 
deviations, 

T is the matrix of cross products of deviations for the total 
sample, 

B is the matrix of between-group cross products of deviations 
of groups from the grand means weighted by group size 
(Cooley and Lohnes, 1962), 


and 

T = B + W. 


*A total of 15 cells occupying the easternmost and southernmost areas was eliminated due to the low fracture trace 
frequency caused by incomplete photo coverage. 
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Figure 29. Rose Diagram Plot of Fracture Trace Patterns for the Kansas 

Study Area 
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The best partition may also be determined by use of the total generalized 
distance, the Mahalanobis D 2 criterion, where D 2 is defined as the sum of 
the distances between multivariate means of all possible pairs of groups, in 
terms of standardized measurements. 

Using principle components, a plot of the eigenvalues of the total correlation 
matrix indicates a gradual decrease in the amount of variation explained by 
each additional component (Figure 31). It was arbitrarily decided to choose 
the 8 component level as the cutoff. A total of 85% of the variation is explained 
by the 8 components in the Kansas data and 82% is explained in the west 
Texas data. 

The two sets of data were processed by the program CLUS, using 2 through 11 
groups. Figure 32 illustrates the plot of the two criteria using the log 
(max | X | / I W | )* algorithm for the Kansas data using 8 components. An 
inflection at the three group level in both criteria is interpreted as significant. 
The six group level also indicates a major inflection of the D 2 * criterion. 

An additional run on the data using six components produced exactly the same 
classification for the 6 group level, but showed a more marked increase in the 
value of both criteria. Figure 33 illustrates the three group classification, 
which in a crude fashion, tends to outline the geology. Group 2 mainly occupies 
the northern part of the area of exposed Illinoian deposits, group three occupies 
the area underlain by the Kansan and Permian deposits and group 1 covers 
areas occupied by a mixture of groups 2 and 3. The six group level (Figure 34) 
contains some isolated members of groups 3 and 5 within the central part of 
the map. These overlie the Coats Anticline, which has a structural closure of 
approximately 250 feet (76 m) and may thus have affected the overlying fracture 
pattern. Examination of the rose diagram patterns in Figure 29 reveals a 
pronounced enhancement of the northeast ray directly over the structure (row 
3, column 3), which may be associated with the northeast trending fault in the 
basement rocks underlying this structure (Cole, 1962). 

Figure 35 illustrates the plot of the two indicator criteria for the west Texas 
data. No pronounced peaks were noted, although a change in slope for both 
criteria occurs at the two and seven group level. The two group classification 
(Figure 36) seems to be related to geologic materials exposed on the surface. 
Gioup 1 tends to overlie areas of Permo-Triassic rocks whereas group 2 
occupies areas of Cretaceous-Tertiary deposits. The 7 group level (Figure 37) 
shows no obvious relation to any of the reef structures. The classification is 
again partly related to lithology; groups 1, 2, and 3 overlie Cretaceous- Tertiary 
deposits, whereas groups 4, 5, and 7 overlie the transition between the two 


‘Based on the grouping of log (max | T | / | W | ). 
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Figure 31. Plot of Eigenvalues versus Principal Components for Both Study Areas 
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KANSAS STUDY AREA 



Figure 33. Results of the 3 Group Classification of Rose Diagrams 
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KANSAS STUDY AREA 



Figure 34. Results of the 6 Group Classification of Rose Diagrams 
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Figure 35. Plot of CLUS Classification Criteria for the Texas Data Using 8 Components, Log 

( |T[ /|W | ) Maximized 
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Figure 36. Results of the 2 Group Classification of Rose Diagrams 
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Figure 37. 


Results of the 7 Group Classification of Rose Diagrams 
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major map units. Group 6 occupies mainly Permo-Triassic rocks. The mis- 
classifi cation of the cells in the southwestern part is probably due to the small 
sample size in this area due to incomplete coverage, producing rose diagrams 
without any preferred rays. 

In summary, classification of the rose diagrams using a multivariate classifi- 
cation scheme produces groupings which are predominantly controlled by 
surface lithologic factors if classification is limited to a small number of groups. 
Active structures (i.e. basement uplift anticlines) may be recognizable because 
their fracture patterns may differ from their immediate neighbors and may be 
isolated by classifications at higher group levels. Passive (reef) "structures 1 ' 
do not create fracture patterns which can readily be isolated from their 
surrounding neighbors by this technique. 

CONCLUSIONS 


Detailed quantitative analysis of fracture trace patterns can be routinely per- 
formed using repetitive techniques and computer algorithms. Cultural 
features can affect the ability to map fracture traces. Fracture trace lengths 
tend to be log-normally distributed. Deviations from log-normality tend to be 
associated with structural closures in both study areas, suggesting that fracture 
pattern may be disturbed over the structures. 

Trend surface analysis may allow extraction of several levels of information 
that may be present in a set of data. Examination of fracture traces by trend 
surface analysis indicates that lithology mainly controlled the frequency and 
log-mean fracture trace length. Frequency was also affected by an operator 
bias, which caused alignment of some of the model contours with flightline 
paths in at least one of the study areas. Higher order surfaces extracted the 
majority of these variations. Residuals in the frequency analysis isolated 
areas of increased fracture frequency in the Kansas area that appeared to be 
associated with either bedrock exposures or with structural culminations. In 
the west Texas area, strong negative residuals appear to be related to reef 
structures. The increase in frequency in the active structure (anticline) and 
the scarcity in the passive structure (reef) suggests either different mechanisms 
for propagation of fractures through these two types of structural discontinuities 
or different stress fields produced above the structures. Analysis of residuals 
for log-mean fracture trace length indicates that in at least one instance, 
fracture traces may be shorter over active structures in the Kansas study area. 
The west Texas residuals map shows some alignment parallel to flightline 
paths; however, a strong negative anomaly (shorter lengths) may be associated 
with a through-going lineament in the area. Both areas show a shortening of 
fracture traces in areas underlain by tectonic structures (anticlines, lineaments), 
possibly due to increased fracturing and subsequent erosion of the fractures. 
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Fracture traces and joints measured in an area underlain by Permian rocks 
(sandstones) coincide in orientation, but there may be large differences in the 
length of the frequency rays. In an area of Cretaceous rocks (limestone), the 
apparent displacement in orientation between joints and fracture traces may 
represent a possible second order shear relationship between the fracture 
trace and jointing directions. 

Analysis of rose diagrams using a multivariate statistical approach shows that 
the basic source of variation is due to differences in surface lithologies, and 
that a lesser amount may be due to deformational effects in an active structure, 
thus changing the fracture pattern. For example, in the Kansas area, an 
anticlinal structure, with a normal fault at depth, was isolated from its surround- 
ing neighbors, whereas in the west Texas area, the predominant effect was 
lithology even in the larger group classifications. 
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FUTURE WORK 


Similar areas should be studied to determine if a valid exploration technique 
has been developed. Additional work also should be carried out to determine 
if lithologic control may be extracted from fracture trace orientations 
summarized as rose diagrams. Conversion to percent rose diagrams may 
achieve this end. 
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APPENDIX A 

SOURCE LISTING OF FORTRAN IV COMPUTER PROGRAM "VECLEN" 
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***************************** VLNOO 305 

VECTOR LENGTH PROGRAM VLN00015 

***************************** VLNOO 0 25 

VLN00035 

THE PROGRAM WAS WRITTEN BY MELVIN POOWYSCCKl OF THE G EDS C I ENCES VLNOO 045 
DEPT., THE PENNSYLVANIA STATE UNIVERSITY, APRIL, 1972 -OR THE VLN00055 
IBM 260/67 COMPUTER. AND WAS MODIFIED IN APRIL, 1974, FOR USE VLN00065 
ON OTHER COMPUTERS HAVING THE EQUIVALENT OF 120K BYTES S TOR AGE . VLNOO 0 75 

VLNOO 085 

PROGRAM SUMMARIZES VECTOR CATA AS FREQ UEN CY-LENGTH D I STR I QUTI ONS . VLNOO 0 95 
UPERATCR SPECIFICS CLASS LENGTH, SIZE OF MAXIMUM CLASS AND DFT ER— VLNOO 105 
MINES THE CIMENSICNS OF THE AREA (GRID CELL) IN WHICH THE DATA VLN001I5 

ARE SUMMARIZED. A VECTOR IS COUNTEO IN A GRID CELL IF ITS MID- VLN0012S 

POINT FALLS IN THE CELL. NG CONSIDERATION IS GIVEN TO VECTOR AZI-VLN0G135 

MUTH. PROVISION IS MADE FOR A LOG BASE 2 TRANSFORMATION IF DESI- VLN00I45 

RED TO ATTEMPT NO K M AL I Z AT I ON OF THE DATA. A TEST FOR NORMALITY ISVLN00155 
MADE BY COMPARING A THEORETICAL DISTRIBUTION USING THE CALCULAT COVLNOO 1 65 
ME A F AND STANDARD DEVIATION OF THE OOSERVED POPULATION AGAINST VLNOO 1 75 
THE DISTRIBUTION CF THE OBSERVED POPULATION UTILIZING THE CHI VLNOO 1 85 

SQUARE CRITERION. EXAMPLES OF PROGRAM OUTPUT ANO APPLICATIONS AREVLNQO 195 
GIVEN IN. POOWYSOCKI, M.h.. 1974, "ANALYSIS OF FRACTURE TRACE VLN00205 

PATTERNS IN AREAS OF FLAT-LYING SEDIMENTARY ROCKS FOR THE OE TEC— VL N00 2 1 5 
TIGS OF EURI ED GEOLOGIC STRUCTURE"; NASA - GODDARD SPACE FLIGHT VLNC0225 
CENTER DOCUMENT X- 52 3- 74-2 C C . DATA ARE P 6 AD FROM CARDS GENERATED VLN00235 
BY VECTOR TRANSFORM PROGRAM (SEE PODWYSOCKI, M.H,, 1974, "FORTR ANVLNOO 24 5 

IV PROGRAMS FOR SUMMARIZATION AND ANALYSIS OF FRACTURE T9 A CE ANO VLN0Q255 
LINEAMENT PATTERNS"; NASA - GSF C DOCUMENT X-644-74-3). VLN00265 

CONTRUL £. TITLE CAROS ARE READ FRLM CARD READER WHILE DATA CARDS VLN00275 
MAY EE READ FROM ANY UN I T DECLARED BY MTAPE2* ON CONTROL CARO 4.VLN00285 


VLNO0295 

ALL NUMERIC INPUT CATA IS RIGHT JUSTIFIED; "I" INDICATES INTEGER VLN00305 
FORMAT, "F" INCICATES FL CATI NG PCINT FORMAT, "A" INDICATES CHA- VLNC0315 
RACTER FORMAT, PROCEEDING NUMBERS INDICATES COLUMNS USED FOR VLN00325 

EACH PARAMETER. TO SPECIFY NCNUSE OF AN OPTION, PUNCH 0 VLN00335 

VLNOO 345 

********CCNTROL CARDS 1 THRU 3 TITLE CARDS VLN003S5 

TITLE WILL BE PRINTED AT THE TUP Of EACH GRID CELL SUMMARIZED VLNOO 365 

( 20A4, #1 -80 I . note: 3 CARDS must BE LSEO; IF ALL 3 ARE NOT USEDVLN0G375 
, BLANK CAROS MUST HE INSERTED IN THEIR PLACE. VLN003Q5 

********CCNTFCL CARO 4 OPTIONS CAPO VLN00395 

X INC= I NCREMENT OF X-AXIS TRAVERSE IN MM, (14,01-4) VLN0040S 

Y I NC=INC«EMFNT OF Y-AXIS TRAVERSE IN MM. ((4,05-8) VLN00415 

XST AFT= ST ART ING POINT FOR X-AXIS TRAVERSE IN MM. (14, *9-12) VLN00425 

YSTART^STARTING POINT FOR Y-AXIS TRAVERSE IN MM. (14.013-16) VLN0043S 

XSTOF=END OF X-AXIS TRAVERSE IN MM. (14,017-20) VLN0Q445 

YST0P-END op Y-AXIS TRAVERSE IN MM. (14,021-24) VLN00455 

NOTE: PROGRAM SUCCESSIVELY SCANS DATA IN MAP GRID CELLS « XCELL • VLNOO 465 
BY 1 YC EL L 1 IN SIZE, INCREMENTING BY * X I NC * UNTIL 'XMAX* > VLN00475 

* X5 T OP * , WHEN • Y I NC p IS INCREMENTED. PROGRAM TERMINATES WHEN VLN00485 

»YMAX» > 1 YSTOP * . NONE OF THE ABOVE 6 VALUES CAN SE NEGATIVE. VLN00495 

AMPS CL ~ MAP SCALE ENTERED AS MILES/MM. ( F 5 ♦ 4.0 25-2 9) VLN00505 

NOTE: WHEN VECTORS are MEASURED ON a l:2400C SCALE MAP AND 0UT-VLN00515 
PUT IS DESIRED IN M IL £ S , * A MPSCL * = • C I 49 (I.E* 1 MM. =.0149 M I LES ) VLNOO 5 25 

OH I NC-NUMER IC AL VALUE OF EACH «X* INCREMENT OF FREQUENCY -LENGTH VLN00535 
HISTOGRAM < I . t . EACH »X»= 2 VECTORS) ( F 5* 2 » 430- 34 ) VLN00545 

SC I NC = FREQUENCY CLASS INTERVAL! DEPENDENT CN DATA TREATMENT VLN00555 

(SEE ♦ NT RAN* BELOW ) « < F5 .2 , 0 35-39 ) VLN00565 

SCLM AX “UPPER CLASS LIMIT OF LAST FREQUENCY-LENGTH CLASS (SEE VLN00575 

• KT R AN 1 BELOW). ( F 5 . 2 • 0 4 C - 44 ) VLN00585 

NH 1ST — PUNCH I FOR FP ECUENC Y-LEN GTH HISTOGRAM IN PRINTED OUTPUT VLN00S95 
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< 11,445) VLN00605 

NPUNCH--PUNCH 1 IF FR EQUENC Y-LENGTH DISTRIBUTION IS OESIRED VLN00615 

ON CARDS < I 1 * #4 £ > VLN0062S 

NST AT PUNCH 1 IF FREQUENCY MCMENTS (l#E. ME AN, STANDARD D E V I AT I ONVLNOO 6 3 5 

, SKEWNESS# ETC • I ARE OE SIRED ON CARDS (11**47) VLN00645 

XC ELL-CELL SI2E (IN MMi> IN X DIRECTION <14, *48-51) VLN006S5 

YC£LL=C£LL SIZE (IN MM.) IN Y DIRECTION <14. *52-55) VLN00665 

NFOLD--PUNCH 1 TO FOLD TAILS OF FREQUENCY-LENGTH DISTRIBUTION TO VLN00675 
SATISFY REQUIREMENTS FOR CHI SQUARE TEST# TAILS ARE FOLDED WHENVLN00685 
EXPECTEO FREQUENCY OF A CLASS 1$ < 0 • 9S (II #456) VLN00695 

IT AP£2=LQG I C AL UNIT FOR READING DATA CARDS GENERATED BY "TRANS- VLN00705 
FORM" PROGRAM (12,457-58) VLN00715 

NT RAN — SELECTS DATA TREATMENT (12*459-60) VLN00725 

PUNCH 1 IF DATA IS TO EE TREATED IN A LINEAR FASHION U#E# VLN00735 

INTERVAL IN MILES AS SPECIFIED IN 'AMPSCLM . NOTE; * NFOLD* C ANVLNOO 74S 

NOT BE 1 IF THIS OPTION IS CHOSEN VLN00755 

PUNCH -1 IF CATA IS TO 0E CONVERTED TO LOG BASE 2 BY THE FOR- VLN00765 
MULAS Z * I 1/LOG1 0(2 > >*L0G1Q< X)+6* WHERE X = VECTOR LENGTH VLN00775 

IN MILES# NOTE: FOLD OPTION MAY BE USED. VLN00785 

1 SCL MAX • AND *SClNC» ARE GOVERNED BY • N TRA N * • FREQUEN CY -LENGTH VLN00795 
CLASSES BEGIN WITH A MINIMUM VALUE OF 0 AND INCREMENT BY « SC I NCVLNOO 80 5 


* UNTIL 1 SCL MAX * IS REACHED. IF A LINEAR SCALE IS USED* « SCLMAXVLNOO 0 1 5 

• IS > THE LARGEST VECTOR LENGTH IN MILES* IF DATA IS TRANSF OR - VLN00825 

MED TO LOGARITHMS#. «SCLMAX‘ IS DETERMINED BY THE CONVERSION OFVLN00835 
THE LARGEST VECTOR LENGTH BY THE ABOVE FORMULA, *SCINC* MUST BEVLN00845 
CHOSEN APPRQPRI ATELY FOR EACH CASE. VLN0085S 

44 44444WCAT A CAROS VLN00865 

VECTOR CATA INPUT FROM VECTOR TRANSFORM PROGRAM VLN00875 

VLN00Q85 

DIMENSION TITLE! 60) *VECLEN( 20C0 ) . XM I D ( 2 000 ) #YMID( 20 00 ) ♦ Z ( 2 000 ) #SCM VLNOO 895 
1 I N ( 4 0 ) » SCM AX ( 401# FNUM <4Q).ZI(40>*AREA<40),DIFF(40) ,FR~ QE X< 40 > * CH LSVLNOO 90S 
2Q< 40 1.0(40) ,FD< 40 ) # FD 2 ( 4 0 1 * FD3( 40 >. FDA ( 40 > *FD5(40) VLN00915 

DIMENSION FFRQ X (40)# FFNUM< 40 > * FCHI SQ ( 4 0 ) VLN00925 

CATA IHX/1HX/,A C/3.321 93 /.l RE AD/5/# I PR l NT/ 6/* I PUNCH/7 / * A 1 / • 099 7926VLN00 93S 
I 8D0 / # A 2/ #04 432 0140 0/# A3/ #00969920 00/* A4/-. 000 09 9620 0/ . A5 /• 0 00 58 1 55VLN 00 945 
200/ VLN0095S 

INTEGER X I NC * Y I NC » XST Afi T * YST A RT * Y MI N * X M I N # XST OP ♦ VSTOP * XM AX . Y MAX , XC VLNOO 96 5 


1 ELL , YCE LL * BOMB 


READ CONTROL CARDS 


VLNOO 9 75 
VLN00985 
VLN0099S 
VLN01005 
VLN01015 
VLN01025 
VLN01035 
VL NO 1 045 


aoMa^o 

RE AC I [READ, 5) ( TITLE! L > #L= 1,60) 

READ (IREAD# 10) XlNC# V I NC , XST ART . Y ST ART , XSTOP , YS TOP * AMPSC L * DH I NC * 
l SC I NC#S CLMA X#NH 1ST, NPUNCH *NST AT , XCELL * YCELL# NFOLD * |T APE2 # MTR AN 
IF < X 1NC.LT #0*OR . Y I NC*LT • 0 #0R# XST ART #LT* 0 *OR. YSTART#LT*0# OR * XSTOP# VLNOI OSS 
LLT.O *CR. YSTOP.LT . O.OR .XCELL #L T . 0 . OR . YCELL# LT * 0 . OR . NFDlD* EQ . 1 . AND VLN01 065 
2.NTR AN. GT # 0 > BOMB— I 
U IFCBCMB) 16.18, 13 

13 WRITE! I PR I NT* I 5> XINC . VINC. XST AR T , YST ART* X STOP ♦ YSTOP ♦ XCE LL . YCELL 
1 .NFOLD# NTfiAN 
GO TO 800 


READ DATA CARDS FRCM LOGICAL UNIT *ITAPE2» 


V L NO 10 75 
VLN0I085 
VLNOI 095 
VLNOi 105 
VLNOI 1 15 
VLNOI 125 
VLNOI 135 


16 DO 25 1-1.5000 

READ (IT APE 2 *20 ♦ END = 30 ) V ECLEN ( l ) , XM ID ( I ) * Y M I D ( I I 
IF! 1^2000 ) 25# 25*22 

22 WR ITE ( 1 PR I NT* 24 > 

GO TC 000 
25 NUM *1 


VLNOI 145 
VLNOI 155 
VLNOI 165 
VLNOI 175 
VLNOI IBS 
VLNOI 195 
VLNOI 205 
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o u L> 


30 

32 

35 

36 

39 

40 

50 

60 

70 


80 


I 


too 


105 

1 10 

120 

130 

140 


DO 35 M=1,NUM 
ZC M|=V£CLENCM> 4 AMPSCL 
IF ( NTRAN ) 32*13*35 

2 C M ) 3 ( ALOGl OCZCM ) )*AC>+6« 

CONT INUE 

NCL ASS= SCLM AX/SC INC+O *5 

IF CNCL ASS-40) 39,39*35 

WRITE! I PR I NT. 38) S CLM AX , SC I NC 

GC TC 800 

DO 70 N=t* NCL ASS 

IFlN-i) 40,40.50 

SCMIM N > = 0 , 

GO TC 6 C 

SCMIKCN )=SCMAXC N— 1 ) 

SCM A X ( N 1 3= SC M I N I N ) + SC I NC 
CONT INUE 

SCAN AND SUMMARIZE VECTORS IN EACH GRID CELL 

DO 700 YMIN-YST ART. YSTCP, Y1NC 
IFf YMIN.GE.YSTOP) GO TC 800 
DO 700 XM1N=XST ART* XSTOP * X INC 
IFtXM|N*GE.XSTOP) <30 TC 700 
DO 80 L — 1 • NCL AS S 
AftEA(LI=0* 

CHISQIL >=0. 

OCL MO* 

OIFF CL I =0 • 

FD(L 1=0 . 

FD2CLM0. 

FD3CL ) = 0 • 

F D4 ( L ) = 0 • 

FD5( L ) = 0 • 

FCHISQt L) =0 • 

FFROXCL 1=0* 

FFNUML 1=0. 

FN'J M (L 1 =0 « 

FREQEXC L ) =0 • 

ZI <L 1 = 0 . 

CONTINUE 
TFNU M=0 • 

XMAXTXMlNfXCELL 
VMAX*YM IN+ Y CELL 
DO 140 1=1* NU M 

1FCXMIDC I 1.GE.XMIN* AND »XMlD(t)«LT . XMA X* AND • YM IO< l > , GC » YM l M*ANO 
DC l I *LT *YMAX> GO TO ICO 
GO TC 140 
I A C= ( 

DO 120 J= 1 * NCL ASS 

IFCZCl AC1 . GE.SCMINC1 ) .AND.ZC ( AC 1 * LT » SCMAXC NCL A SS 1 1 GO TD 1 05 
WRITEC I PR INT.102) I AC .ZC I AC ) ♦ SCM INC 1 ) • SCMAXC NCL ASS I 
GO TO 8 CO 

IFC Z C 1AC1 .GE. SCM IN ( J1 • AND.ZC I AC ) .LT.SCMAXC J) 1 GO TO 110 

GO TO 120 

NT YPE=J 

GO TC 1 30 

CONTINUE 

FNUMCNT YPE1 =FNUM< NTYPE ) + 1 , 

CONT INUE 

DO 150 N= 1 * NCL ASS 
TFNU> = TFNUM + FNUH(N I 


VLN01 215 
VLN01 225 
VLN01 235 
VLNQ1 245 
VLN01255 
VLN01 265 
VLNOl 2 75 
VLNOl 285 
VLN01295 
VLN0130S 
VLN013I5 
VLN0I325 
VLN01 335 
VLN0134S 
VLN0135S 
VLN01 365 
VLN013T5 
VLN01 365 
VLN01 395 
VLN01405 
VLN01415 
VLNOl 425 
VLN01435 
VLNOl 445 
VLNOl 4 55 
VLNOl 465 
VLNOl 475 
VLN0148S 
VLN01495 
VLNOl 5 05 
VLNOl 515 
V L NO 1525 
VLN0153S 
VLNOl 545 
VLNOl 555 
VLNOl 565 
VLN01575 
VLNOl 585 
VLNOl 595 
VLNOl 605 
VLNOl 615 
VLNOl 625 
VLN01635 
VLNOl 645 
YMl VLNOl 6 55 
VLN01665 
VLNOl 675 
VLNOl 685 
VLNOl 695 
VLNOl 70S 
VLNOl 715 
VLNOl 725 
VLNOl 7 35 
VLNOl 745 
VLN01755 
VLN01765 
VLNOl 775 
VLN01785 
VLNOl 795 
VLNOiaOS 
VLNOl 815 
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150 CONTINUE 

NC~< SCLMAX-Q* > /SC INC+O *5 

ELIMINATION OF LONER EMPTY CLASSES 

DO 161 JK=1,NCLASS 
I F < F MJ M ( JK) ) 161*161* 162 

161 CONTINUE 

162 JKL=JK 

ELIMINATION OF EMPTY UPPER CLASSES 

DO 163 JK=1 ,NCLASS 
KH = ( NCLASS-JK) * I 
IF C FNUM <KM>) 163*163* 164 

163 CONTINUE 

164 JKH^KH 

CALCULATE STATISTICAL MOMENTS FOR EACH GRI C CELL 
MAXCLS- 1 

DO 1 66 M-JKL«JKH 

IF (FNUM (MAXCLS >-FNUM< M 1 I 165* 165* 166 

165 MAXCLS=M 

166 CONTINUE 

CLSMCP= (SCMIN(MAXCLS) *5 CM AX I MAXCLS) 1/2 
DO 170 M$= JKL . JKH 

DIMS ) = ( ( SCM 1 N( MS >+SCMAX( MS) >/Z-CLSMDP >/SCI NC 
170 CONT INUE 
SD-0 i 
SD2 — 0 • 

SD3=0. 

SD4=0. 

SD5= 0* 

DO 175 1= JKL* JKH 

FD( 1 ) — FNUM < l l*OC I > 

SD=S C ♦F C( l I 

FD2( K ) -FNUM ( I ) 4 ( D ( I )44£> 

SD2=S02+FD2< I 1 

FD3< l)=FNUMU )*< D< 1)4*3) 

SD3 = S03 +FD3 ( I ) 

FD4< I)=FNUM(I )4<0( I >4*4) 

SD4- SD4 4F D4 ( I > 

FD5 < I > = FNUM ( I > * < ( D ( I )-l*)+*4J 
S05“SD5*FD5< I ) 

175 CONTINUE 

GCK=SD4-<4**SD2)M e*4SD2)-( 4*4SD) 4-TFNUM 

AMOM l = SD /T FNUM 

AMCM2=SC2/TFNUM 

AMOM 3= S 03/ T FNUM 

AMC M4 — SD4/T FNUM 

TMC M 1 = AMOM1 4SC I NC 

TMOM 2= ( SC INC*42)*( A MOM 2— C (AMO Ml )4*2> ) 

T MOM 3— (SCI NC 44 3 >4< AM DM 3 * ( 3 * * A MO M2 4 A MOM 1 )♦( 2 • 4 ( AMOM 1443) ) ) 
PT4A~< AM0M4-C 4 . 4 A M CM3 * AMO Ml )♦ < 6 • 4 ( AMOMl 442 )*AMOM2> ) 
PT4R«PT4A- ( 3*4 ( AM0M1444 ) ) 

TMQM4=( SCINC444 ) 4 ( PT4B ) 

XBARvCLSMDP +TMOM 1 
VAR=TMCM2 
STDV3SORT ( VAR) 

RT 01 — (TMCM3/( ( SORT (VAR 1)443)1 


VL NO 18 2 5 
VLN01835 
VLN01 845 
VLIM01 855 
VLN01 865 
VLN01875 
VLN01 835 
VLN01 695 
VLN01 905 
VLN019I5 
VLN01925 
VLN01935 
VL NO 1 94 5 
VLN0195S 
VLN01965 
VLN01975 
VLN01985 
VLN0199S 
VLN02005 
VLN020 1 5 
VLN02025 
VLN02035 
VLN0204S 
VLN020S5 
VLN02065 
VLN02075 
VLN02085 
VLN02095 
VLN02105 
VLN021 1 5 
VLN021 25 
VLN021 35 
VLN021 45 
VLN021 55 
VLN021 65 
VLN02175 
VLN021 85 
VLN0219S 
VLN02205 
VLN02215 
VLN02225 
VLN02235 
VLN0224S 
VLN02255 
VLN02265 
VLN02275 
VLN02285 
VLN02295 
VLN02305 
VLN02315 
VLN02325 
VLN0233S 
VLN02345 
VLN02 355 
VLN0236 5 
VLN0237S 
VL NO 2 3 6 5 
VLN02395 
VLN02405 
VLN024 1 5 
VLN02425 
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200 


210 

220 

230 


24 0 

250 

260 

2 70 
280 


290 


C 

c 

c 

300 


302 

305 


310 


B2 »<TF0M4/<TMCM2**2 ) > 

CALCULATE CHI SQUARE CONTRIBUTION 0 EXPECTED FREQUENCIES FOR EACH 
CLASS ASSUMING A NORMAL DISTRIBUTION WITH A XBAR & STOV OF THE 
OBSERVED POPULATION 

NDF=C JKH-JKL)-2 

DO 200 J— JKL* JK H 

Zl < J ) = < SCM AX C J ) — xe AR ) /STD V 

CO NT INUE 

TO I FF=0 . 

TFRQEX=0. 

TCHISQ^O* 

I JK = 0 
IKC= JKL-l 
AREAUKO-O - 

00 290 1 = JKL * JKH 

1 I V= l-l 

IF ( Z I < I )) 21 0*220* 220 
AVT- — Z I ( I I 
GO TO 230 
AVT-2IC I 1 

AREA i I )=i*+AVT4< A1 + A V T* ( A 2 + A VT* { A 3 + A VT * t A4 +A V T* A5 ) • ) > 

AREAII )=0#5/<AREAt I )+*B) 

IF <Z K l n 240*250.250 
OIFF < I |= ARE MI I — AR EA{ I IV) 

GO TO 280 
I JK= IJK + t 

IF<IJK-1) 26 0*260*270 

DIFF< I >=1*-AREA< I ) — AR E A ( I IV) 

GO TC 2eo 

DIFF I I ) = AREA( I IV) — ARE A f I) 

TD IF F“T CIFF+D1FFI I ) 

FREQEX < I )-= D IFF ( I ) ATPNUM 
TFR 0EX= TFRQEX+ FREQ £ X< I ) 

CHISQI I )=< ( FNUM( l ) — F REO EX ( I ) )**2)/FREQEX(I) 

TCHI SO=TCHI SQ + CH ISQ( I) 

CONT INUE 

CHIPR0-PR8CHI ( TCHISQ* NDF) 

NFL A G=0 
MFLAG-0 

IF I NFQL 0)495*495* 300 

FOLD LOWER TAIL OF DISTRIBUTION IF REQUIRED 

DO 302 KP= JKL . JKH 
FCH 1 SQ ( MP ) sCH ISO C MP ) 

FFRQK(MF) =FREQ E X ( MP ) 

FFNUMIMP) =FNUM ( MP ) 

CONTINUE 

00 310 LL- JKL • M AXCLS 
IF<FFRQX<LL )-0*9S> 3O5*31Ct310 

FFRQ X (LL+ 1 ) =FFRQX ( LL ) ♦ FFRQX < LL+ 1 ) 

FF NUM< LL+ 1 >=FFNUM(LL> +FFNUMI LL + 1 ) 

FCH I SQILL )=0. 

JKL 1*LL+1 
NFLAG^i 
CONT INUE 

JKQ= < JKH-MAXCLS 1+1 

FOLD UPPER TAIL OF DISTRIBUTION IF REQUIRED 


VLN02435 
VLN02445 
VLN024SS 
VLN02465 
VLN02475 
VLN0243S 
VLN0249S 
VLN02505 
VLN025 I 5 
VLN02525 
VLN02S3S 
VLN02545 
VLN02SSS 
VLN02S6S 
VLN0257S 
VLN02585 
VLN02595 
VLN0260S 
VLN0261 5 
VLN02625 
VLN02G35 
VLN0264S 
VLN02655 
VLN02665 
VLN0267S 
VLN0268S 
VLN02695 
VLN0270S 
VLN0271S 
VLN02725 
VLN02735 
VLN02745 
VLN02755 
VLN02765 
VLN02775 
VLN0278S 
VLN02795 
VLN0280S 
VLN028I5 
VLN02Q25 
VLN02835 
VLN02845 
VLN02855 
VLN0266S 
VLN02375 
VLN02685 
VLN02A95 
VLN0290S 
VLN02915 
VLN02925 
VLN02935 
VLN02945 
VLN02955 
VLN0296S 
VLN02975 
VLN02985 
VLN02995 
VLN03005 
VLN03015 
VLN03 0 2S 
VLN03035 
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riftn 


OF THE 


OO 320 LH=1.JKQ 
KHH= (JKH-LHJ+l 

IF (FFRQX(KHH)“0*95 I 315*320*320 
315 FFRQ> (KHH- 1 ) = FFRQX (KHH-1 I 4-FFRQX < KHH ) 

FFNUM (KHH— 1 >=FFNUM (KHH- 1) ♦FFNUM( KHH ) 

FCH I SO I KHH ) =0 « 

JKHI sKHH- 1 
MFLAG=1 
320 CONTINUE 

IF < NFL AG ) 325,325*330 
325 J2=JKL 

(30 TO 335 
330 J2*JKL1 

335 IF ( MFL AG) 340*340,345 
340 J3-JKH 

GO TO 350 
345 J3=JKHl 
350 TFFR QX=* 0 ♦ 

TFCHSO^C. 

OO 355 J|sJ2,J3 
TFFRQ#wTFF«QX4FFRqX< J I I 

FCH ISO < J< )=< I FFNUM I JI >~FFRQX< JI ) ) 442 > /FFRQ X< J I ) 
TFCHSQ^TFCHSQtFCHlSQC JI ) 

355 CONTINUE 

NFDFM J2-J21-2 

FCHPFB=PR 0CHICTFCHSQ* NFDF ) 

OUTPUT 

495 WR I T E ( I PR I N T , 50 0 ) ! T I TLE < L } *L* I , 6 0 > 

NRQW=< YMIN4YINC INC 
NCOL=*(XMIN*XINC)/XINC 

WRITE I I PR I NT* 5101 NRO * *NCOL • XMI N* XMA X * YM I N * YM AX 
WRITE! I PR I NT , 520 ) 

NXERIfeO 

DO 630 I J — J Klw * JKH 

WRITE! I PR I NT *530) I J,SCMIN( I J |,SCMAX( I J I »FR£QEX( ( J) ,CWI$Q( I J) * 
lFNUM(U) 

IF INFOLD *EQ.l * AND . NFL AG .EO . 1 .AND, J2.EQ.I J) WRITE! I PRINT, 53 21 
1FFROXI Ivll * FCHI S0< I J),FFNUM< I J ) 

IFINFOLO.EQ *1 , AND . MFL AG ,EO ♦ 1 , AND, J3.EQ* I J) WR I TE I I PR I NT » 53 2> 
IFFRQXI IJ> ,FCM1S0< I J ) * FFNUMI IJ) 

IF INPUN CH-1 I 560*540*560 

54 0 WRITE! I PUNCH, 550) NROW ,NCOL. I J,SCHIN( i J> , SC MAX! I J) , FNUM { I JI 
560 IF (NHIST-1) 630*570,630 
570 NUMXic FNUM ! U) /OH INC 
IF(NUMX) 600*550*600 
560 WRITE! I PR I NT , 590 ) 

GO TC 630 

6C 0 IFINUMX-70) 620.620*610 
610 NUMX*70 
NX ERR- I 

620 WRITE! IPRINT.625I I IHX * I UKA = 1 * NUM X I 
630 CONTINUE 

WRITE! I PR I NT #640) TFR Q£X* TCHI SO* TFNUM # OH I NC 

IF INFOLC* EQ . 1 *AND*NFLAG #EQ* 1 .OR .MFL AG .EQ# 1 > WRITE !IPRINT,64S) 
1TFCHSQ 

IF! NXERR.EQ.U WR ITE< IPR I NT . 6 50 ) 

WR !T £ I I PR l NT *65 2 ) NOF,CH|PR0 

IF < NF OLO #EQ « I • ANO* NFL AG »£G • 1 • OR *MFL AG *EQ« 1 ) WRITE (IPRINT.654) 



BEPRODUCIBILnY 
ORIGINAL PAGE IS 


VLN03045 
VLN03055 
VLN03065 
VLN03075 
VLN03085 
VLN03095 
VLN03105 
VLN031 1 5 
VLN03125 
VLN03135 
VLN03145 
VLN031 55 
VLN03165 
VLN03175 
VLN03I8S 
VLN03 i 95 
VLN0320S 
VLN03215 
VLN03225 
VLN03235 
VLN03245 
VLN03255 
VLN03265 
VLN03275 
VLN03295 
VLN0329S 
VLN03305 
VLN03316 
VLN03325 
VLN03335 
VLN03345 
VLN03355 
VLN03365 
VLN0337S 
VLN03385 
VLN03395 
VLN03405 
VLN03415 
VLN03425 
VLN03435 
VLN03445 
VLN0345S 
VLN03465 
VLN03475 
VL NO 3 4 Q 5 
VLN0349S 
VLN03505 
VLN03515 
VLN0J525 
VLN03535 
VLN03545 
VLN03555 
VLN03565 
VLN03575 
VLN03S8S 
VLN03595 
VLN03605 
VLN036 IS 
VLN03625 
VLN03635 
VLN0364S 



c 

c 

c 

c 

c 

c 

c 

c 

c 


I NFDF *FCHPRB VLNC3655 

WRITE! L PR 1 NT *655 1 MAXCLS» CL SHOP * F NUM ( H AXCL S I VLN03665 

WRITE ( IPR I NT *660) X8 AR , V AR • S TD V « RTB l *B2 VLN0367S 

IF (NSTAT.EQ.I) WRITE! [PUNCH, 6651 NROW*NCOL , XH IN , XMAX* VMI N* VMAX » VLN03685 

1 XBAR * STDV, RTB1 ,82 VLN03695 

IF f GCK-S05.GT • 10* 1 WR [TEC IPR ! NT ,670 > GCK.SD5 VLN0370S 

700 CONT INUE VLN03715 

5 FORM AT ( 20 A4 ) VLN03725 

10 FORM AT! 6 14 *FS.4« 3F5.2 • 31 l ,2 14 • I 1 » 21 2) VLN03735 

15 FORMAT ( * l * * • ONE OR MORE OF THE FOLLOWING DO NOT CONFORM TO PROGRAM VLN03745 

1 LIMITATIONS**/* * *4X ,* XINC •* 4X. • YINC * *2X* *X$TART« ,2X* *Y START • VLN03755 

2 * 3 X * <XSTOP* *3X, * VST OP* *3X * • XCELL 1 ,3X, • VCELL • • 3X , * NFOLO* * 3X , VLN03765 

3 • NTR AN • */ • •* IOI0*/*O * *30X* • JOB ABORTEO * I VLN03775 

20 FORM AT ( T 25 ♦ F6 • t » T5 6 • 2F7 ♦ t > VLN03785 

24 FORMAT ( • 1 * * *MQRE THAN 2000 VECTORS IN DATA SET • . / 1 0 • » 30X , • JOB ABORVLN03795 
1TED* I VLNO3805 

38 FORMAT { *1 ♦* *NUM0ER OF FRE OUENCV— LENGT H CLASSES EXCEEDS PROGRAM LIMVLN038I5 
, 1 I T OF 40 WHERE: SCLMAX/SC INC—NCLASS '/ • • • *T69*F6.2, • / • *F 5* 2 • • = VLN03825 

2** I5*//*0».30X* • JOB ABORT ED * 1 VLN03835 

102 FORMAT! *1 •* 'VECTOR W**I6,* LENGTH OF**FlI*5** EXCEEDS CLASS L I M I TS VLN0384 5 
1 OF* *2F 12.5 ♦/• 0* *3CX, * JOB ABORTED*) VLN03855 

500 FORMAT! • I ♦ ♦ 3! T 3 0 * 2 CA4 * // > ) VLN03865 

510 FORM AT( *0 * # *ROW * * 14. * * COLUMN • • I 4 ,T 30 * • ( • * 15 *• < X < » ,15 *• ; VLN03875 

1*,I5 , • < V < **I5 ♦ *)*) VLN03885 

52 0 FORMAT! *0* *T9. • LOWER * *2X* *UPPER ' ,T42, *CHI • ./*T9 * * CL ASS* • 2X * *CLAS S • VLN03895 
1 *5X* * EXPECT EO • *7X. * SQUAR E •• 5X *• OBSERVED * ./ *T2, 'CLASS* *2X * • LI MIT « » 2VLNQ3905 
2X , •LIMIT** 5X* • FREQUENCY • *5X , «CONTRIB. » *4X# 'FREQUENCY *. 9X • * 0BSERVE0VLN039 1 5 

3 FREQUENCY HISTOGRAM*,//) VLN03925 

530 FORMAT!* • « 2X, I 2 * 2 X , F5 . 2 * 2X , F 5 • 2 « F 7 . 2 * 8X * F 8. 2 , 9X.F4.0) VLN03935 

532 FORMAT! • + ■ .T27. • ( * . F5 . 2 , * ) • * T4 3 . • ( • *F6. 2. « > •* T56. * ! * ,F3. 0* • > * ) VLN03 94S 

550 FORMAT! 315* 2FI 0.2 ,F 12.2) VLN03955 

590 FORMAT! *4* ,T61 **>• ) VLN03965 

625 FORMAT! •♦* ,T6l *•>* ,70A 1 ) VLN03975 

6 40 FORMAT! • * *T20*7< '-•)«T36*7("-')»TS3« » **/«TlO. 'TOTALS* *T | 9 • F8 • VLN039 85 

12*T35*F€.2.T50,F6.0 ,T70, 'EACH "X** = **F10.2.* VECTORCSM) VLN03995 

645 FORMAT!* * * TJ4 , • ( • , F6 • 2 , * ) * ) VLN04005 

650 FORMAT! *0 • ,T70* * ONE OR MORE FREQUENCY CLASSES EXCEED HISTOGRAM LIMVLN04Q15 

1 IT S * ) VLN0402S 

652 FORMAT! *0 *. ‘CEGREES OF FREEDOM (NON-FOLDED) = **I4,*; CHI SQUARE VLN04035 

1 PRO0 ABILITY = *,E10.4) VLN04045 

654 FORMAT! *0 DEGREES OF FREEDOM (FOLDED) = **I4**; CHI SQUARE PROBVLN04055 

IABILITY = **E10.4> VLN04065 

655 FORMAT! *0* *T5Q* • MODAL ST A T I ST ICS * */T50 • l 6< *-•>»//. T20 *• CLASS ** T40 , VLNQ4 0 7 5 

I • M I 0 FOi NT * * T60, * OB S « FREQUENCY • •// »T21 .1 2.T 40* F 7 . 3 * T64 *F6 .2 ) VLN040 8 5 

660 FORMAT! *0 * * T50 * • ST AT I ST ! C AL MCMENTS ',/,T50*I9!*-«> *//T20 *• AVER AGE • VLN04 09 5 
1 *T40* 'VARIANCE* *T60** STANCARD DEV I AT I ON • *T 65* • ROOT Bl «*TI05.* BVLN04 10 5 

2 2 * »//*T20*F8*3 , T40 *FB • 3* T6S *F B.3 * T85*F8 .3* T 1 05* F8 . 3 ) 

665 FORM AT( 614* 4F6.3 > 

670 FORMAT! *0* , * Gfl AM-CH ARL! ER CHECK = *»F1S.4,« JSUM = *,F15.4) 

800 STOP 
END 


FUNCTION PRBCHI (CHISQ. IDF) 


VLN041 15 
VLN04125 
VLN041 35 
VLN0414S 
VLN0415S 
CHI 0005 
CHI 0015 
CHI 0025 
CHI 0035 
CHI 0045 
CHI 0055 


WRITTEN BY H.D. KNOBLE t F. YATES BORDEN* THE PENNSYLVANIA STATE 
UNIVERSITY* 1966 

THIS FUNCTION COMPUTES BY THE APPROXIMATIONS ON PAGE 941 OF 
••HANDBOOK OF MATHEMATICAL F UNCT l ONS •• * U. S. DEPT. OF COMME RC E * 1 96 4 • CHI 0065 
GIVEN A VALUE OF C H I — SOU ARE AND ITS DEGREES OF FREEDOM* FUNCTION CHI 0075 
PRBCHI COMPUTES THE PROBABILITY OF A GREATER VALUE OF CH I— SQUARE. CHI 0085 
THE t ( ARGUMENT ) FUNCTION IS COMPUTED BY FORMULA 26.2*1, P. 931* CHI 0095 
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TR fproducibility of the 
ORIGINAL PAGE IS POOR 


INTEC ER TEST 

ALL REAL*8 ARGUMENTS CHANGE C TO DOUBLE PRECISION BY M • PD Dtf YSOCK I < 
DOUBLE PRECISION DSQRT .QEXP ,ARG.SCHlSO» XPLEVL 
DOUBLE PRECISION Q • R * S ♦ T . L * V • V9 • PRCB * S2PI * ZOOS • APPROX 
DATA S2PI/2. 5066282000/ 

Q( ARG) = ( DE XP( - ARG*ARG*0 *5 )/ 2. 5066262000)*! T*< O *3193 BX 5D0 0*T* 

1 < -0 * 356 5638 000 *T*( 1 .7E 14 7 8DOO+T * ( -1.821 256000+ ( l • 330 274D 0 0*T > > > ) i 
XPL = 2 • 5762 3596D 0 0 
PRBCHl=0.0 

IF (CHI SO .LT .0.0 ) RETURN 
IF ( IDF.LE.O > RETURN 
100 SCHISQ=CHISQ 
S= 1 . 0 
V=IOF 

V9=2 AO/FLO ATI 9* I OF ) 

U=-SCHl 50*0 . 5 

SCH I S Q= DSQRT (SCHISO) 

IF ( C A0S< U ) .LT .174*6) GO TO 110 

174,6 IS THE LARGEST ARGUMENT THAT EXP WILL TAKE. 

PROBsO. 0 
GO TO 240 

CHECK FOR DEGREES OF FREEDOM GREATER THAN 100 OR GREATER THAN 30 

110 IF < IDF . GT .10 0) GO TO 200 
IF ( I OF »GT . 30 ) GO TO 170 

DEGREES OF FREEDOM LESS THAN OR EQUAL TO 30 

PRoe*o.o 

TEST*M0C<IDF,2I 

IF ( TEST . NE • 0 ) GO TO 140 

EVEN DEGREES OF FREEDOM ** LESS THAN OR EQUAL TO 30 ** FORMULA 

2£.4.S* PAGE 941 

IRANCE-I l OF— 2 1/2 

IF ( IRANGE.EQ.O > GO TO 130 

DO 120 I=1,IRANGE 

IR= I * I 

S=S* IR 

120 PROB-PRCB+SCHI SQ*4 IR/S 
130 PROBED EXP (U ) * ( 1 .0 +PR OH ) 

GO TO 230 


C 

c 

c 

c 


ODD OEGREES OF FREEDOM 
26.4.4, PAGE 941 


** LESS THAN OR EQUAL TO 29 ** -DRMULA 


140 IR ANGE— I IOF-1 >/2 

IF (IRANGE.EQ.O) GO TO 16C 
DO 150 1 = 1* IR ANGE 

IR= I +1-1 
S=S*IR 

150 PRQB«PRC04SCHISQ** IR/S 

160 T=1.0/( 1 .0*0.2 3164 19D0C4SCH I SQ) 

PROB^2.0*(Q(SCHISQ>)*2 .0 * I DEXP ( U ) /S2P I >*PROB 
GO TO 2 30 


CHI 

0 105 

CHI 

0115 

CHI 

0 125 

CHI 

0 135 

CHI 

0 145 

CHI 

0155 

CHI 

0165 

CHI 

0175 

CHI 

0185 

CHI 

0195 

CHI 

0205 

CHI 

0215 

CHI 

0225 

CHI 

0235 

CHI 

0 245 

CHI 

0255 

CHI 

0265 

CHI 

0275 

CHI 

0285 

CHI 

0295 

CHI 

0305 

CHI 

0315 

CHI 

0325 

CHI 

0335 

CHI 

0345 

CHI 

0355 

CHI 

0 365 

CHI 

0375 

CHI 

0385 

CHI 

0395 

CHI 

0405 

CHI 

0415 

CHI 

0425 

CHI 

0435 

CHI 

0445 

CHI 

0455 

CHI 

0465 

CHI 

0475 

CHI 

0405 

CHI 

0495 

CHI 

0505 

CHI 

0515 

CHI 

0525 

CHI 

0535 

CHI 

0545 

CHI 

0555 

CHI 

0565 

CHI 

0575 

CHI 

0585 

CHI 

0595 

CHI 

0605 

CHI 

0615 

CHI 

0625 

CHI 

0 6 35 

CHI 

0645 

CHI 

0655 

CHI 

0665 

CHI 

0675 

CHI 

0685 

CHI 

0695 

CHI 

0705 
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C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


CHI 

********** GREATER THAN 30 DEGREES OF FREEDOM ********** CH j 

AN APPROXIMATE VALUE OF CHISQ IS FIRST COMPUTED THEN COMPARED WltH CHI 

the given chi so* if the approx* value is greater Than the given chi 

VALUE* 0< CH ISO* 1 OF I IS RETURNED AS *996. CHf 

*********************************** *********************** ********** CHI 

FOR GREATER THAN 30 AND LESS THAN OR EQUAL TO 100 DEGREES O- FREEOdMCHl 
THE APPROX. VALUE OF CHISQ AT THE .90S LEVEL IS COMPUTED BY FORMULA CHI 
26.4.17. PAGE 941. THE SIGN OF X(P> IN THE FORMULA WAS CHANGED 
FROM * TO - TO ALLOW COMPUTATION OF CHISQ AT THE .995 LEVEL RATh£R 
THAN THE .005 LEVEL AS IS THE CASE WHEN THE SIGN IS ♦ . 


170 


180 

190 


APROX=< ( 1 ,C-V9-XPL*DSCRT< V9 ) >**3>*V 
IF C APROX.LE.CH ISQ ) GO TO 180 
GO TC 210 

V=C (CHI SQ7V )**0 .33 233 333000- ( I.0-V9) >/OSQRT( V9> 
T= I .0/1 1 .0*0.23164 19D0CAV ) 

PROe*Q( V> 

GO TO 2 30 


CH I 
CHI 
CHI 
CHI 
CHI 
CHI 
CHI 
CHI 
Chi 
CHI 

Chi 


GREATER THAN 100 DEGREES OF FREEOOM. THE APPROX. VALUE OF CHISQ 
IS COM PUT EO BY FORMULA 26.4.16, PAGE 941. THE SIGN OF X(^) WAS 
CHANGED FOR THE SAME REASON AS ABOVE, 

2 00 APROX=( (*“ XPL4DSQRT (V + V-l.O) )**2)*0.5 

IF C APRCX.LE. CHISQ ) GO TO 220 
210 PR08**0-99S 
GO TO 2 40 

220 V-DSQRT 12 »0 00*C H ISO— DSORT ( 2. 0*V-1.0) 

GO TO 190 

230 IF (PROE.GT .0.9S5I GO TO £10 
240 PRBCHI=PROB 
RETURN 
END 


CHI 

CHI 

CHI 

CHI 

CHI 

CHI 

CHJ 

CHI 

CHI 

CHI 

CHI 

CHI 

CHI 

CHI 

CHI 

CHI 


0715 
0726 
0735 
0746 
0755 
0 765 
0775 
0785 
0795 
0805 
0815 
0625 
0835 
0845 
0655 
0665 
0875 
0885 
0895 
0905 
0915 
0925 
0935 
0945 
0955 
0965 
0975 
0985 
0995 
1005 
1015 
1025 
1035 
1045 
1065 
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